Functional Integral for Ultracold Fermionic Atoms 
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We develop a functional integral formalism for ultracold gases of fermionic atoms. It describes 
the BEC - BCS crossover and involves both atom and molecule fields. Beyond mean field theory 
we include the fluctuations of the molecule field by the solution of gap equations. In the BEC 
limit, we find that the low temperature behavior is described by a Bogoliubov theory for bosons. 
For a narrow Feshbach resonance these bosons can be associated with microscopic molecules. In 
contrast, for a broad resonance the interaction between the atoms is approximately pointlike and 
microscopic molecules are irrelevant. The bosons represent now correlated atom pairs or composite 
"dressed molecules". The low temperature results agree with quantum Monte Carlo simulations. 
Our formalism can treat with general inhomogeneous situations in a trap. For not too strong 
inhomogeneities the detailed properties of the trap are not needed for the computation of the 
fiuctuation effects - they enter only in the solutions of the field equations. 

PACS numbers: 03.75.Ss; 05.30.Fk 



I. INTRODUCTION 

Ultracold fermionic atoms can exhibit both the phe- 
nomena of a Bose-Einstein condensate (BEC) P, Q of 
molecules and the condensation of correlated atom pairs 
similar to BCS-superconductivity [1, Recent experi- 
mental progress [1^ @, 0) Q in the crossover region reveals 
the universality of the condensation phe nomenon, as an- 
ticipated theoretically [1, [13, HH [12, 

In this paper we develop a systematic functional in- 
tegral formulation for the treatment of the equilibrium 
state ^ of ultracold fermionic atoms. We discuss in detail 
how to arrive at a formulation that treats the fermionic 
fluctuations of unbound atoms and the bosonic fluctua- 
tions of the molecule or di-atom field on equal footing. 
An approach based on a Hubbard- Stratonovich trans- 
formation is an ideal starting point for a unified inclu- 
sion of fiuctuations of molecules or "Cooper pairs". We 
show how this formalism can be implemented in practice 
in a self-consistent approximation scheme. We carefully 
discuss the renormalization procedure that is needed to 
absorb the ultraviolet divergences in this nonrelativistic 
quantum field theory. The result is an effective low en- 
ergy formulation which is insensitive to the microphysical 
cutoff scale A. 

We concentrate on dimensionless quantities by mea- 
suring all quantities in units of the Fermi momentum kp 
or the Fermi energy ep = kp/2M. The inverse Fermi 
momentum is the most important length scale in the 
problem, related to the total number density of atoms 
by kp — (STT^n)^/'^. It measures the typical interparticle 
spacing. We show that the atom density does not appear 
as an independent parameter in the computations which 
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can be performed in terms of dimensionless ratios. This 
renders our formalism highly universal, since the results 
of experiments with different atoms and densities can be 
related by simple scaling laws. 

For this purpose, we introduce three dimensionless pa- 
rameters that characterize the crossover problem effi- 
ciently: First, the "concentration" c describes the ratio 
between the in-medium scattering length and the aver- 
age distance between two unbound atoms or molecules. 
Its inverse, c~^, smoothly connects the weakly coupling 
BCS regime (c < 0, |c| ^ 1) with the BEC regime 
(c > 0, |c| <JC 1). The perhaps most interesting region 
is the "crossover regime" in between, \c~^\ < 1. Sec- 
ond, a Yukawa or Feshbach coupling /i^ characterizes the 
interaction between atoms and molecules. The third pa- 
rameter is the temperature in units of the Fermi energy, 
T = T/ep. No further details of the "microphysics" are 
needed for the macroscopic quantities. In this sense the 
description becomes universal. In the language of quan- 
tum field theory or critical phenomena the parameters c 
and describe relevant couplings for the long distance 
physics. 

Our formalism is well suited to describe all regimes of 
coupling and temperature, including the superfluid phase 
with broken symmetry. As a first application we compute 
the phase diagram for the crossover problem as described 
by the dependence of the critical temperature on c. The 
function T(c^^) only depends on the value of hcj, and 
shows universal limits for large |c~^|, i.e. in the BCS 
or BEC regime, respectively. In the crossover region the 
dependence on the Yukawa coupling is strongest. Nev- 
ertheless, we find that for the two limits ^ Q (narrow 
resonance) and oo (broad resonance) the crossover 

becomes independent of h^. For broad Feshbach reso- 
nances the Yukawa coupling becomes an "irrelevant" pa- 
rameter. 

In the narrow resonance limit for ~* the bosonic 
degrees of freedom can be associated with microscopic 
molecules. In this limit the molecule fluctuations can 
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be neglected and mean field theory becomes valid. In 
contrast, for ^ oo the microscopic molecule degrees 
of freedom play no role and the model is equivalent 
to a purely fermionic model with pointlike interactions. 
Nevertheless, effective molecular bound states ("dressed 
molecules") become a crucial feature of the crossover 
physics. 

For an actual comparison with experimental results 
one needs to relate the "universal parameters" to ex- 
perimental observables, in particular the strength of the 
magnetic field. This is done in 15], where we connect 
the scattering properties of atoms with the values of our 
universal parameters. The present paper may therefore 
be viewed as the theoretical basis for the more detailed 
comparison with experiment in [isj . Here we concentrate 
on the conceptual and formal developments. 

A further aspect of universality concerns the geometry 
of the trap. We will present a systematic "derivative 
expansion" which computes the effects of fluctuations 
and interactions independently of the shape of the trap. 
The details of the trap only enter at the end through 
the solution of effective field equations in presence of a 
trap potential. The systematic character of the deriva- 
tive expansion allows for a quantitative estimate of the 
reliability of simple approximations, like the "local den- 
sity" resp. "Thomas- Fermi approximation" or the "local 
condensate approximation" frequently used in literature 
[l6, 17, 18, 19, 2Q\. The effective field equation for the 
condensate has the same status as the (time independent) 
Gross-Pitaevskii equation [2l|, [HI for a Bose-Einstein 
condensate. 

Our approach basically relies on two ingredients: The 
functional integral and the Hubbard- Stratonovich trans- 
formation or partial bosonization. The partial bosoniza- 
tion permits us to formulate the problem microscop- 
ically as a Yukawa theory, thereby allowing to deal 
with nonlocal interactions. This route is also taken in 

[n, [13, n, m, m, Hi m, ii. Hi [m, m, m . The power of 

the functional integral techniques, however, is so far only 
marginally used. A first attempt in this direction was 
made by Randeria [O, [13]. Later, other approaches em- 
ployed this concept more as an argumentative tool than 
as a method for concrete calculations [1^, [s^, [U, [35[ ^ . 

Beyond the systematic functional integral formulation 
and the emphasis on the universal aspects of the phase 
transition our work extends previous results by the sys- 
tematic inclusion of the molecule fluctuations. These 
fluctuations are important for the quantitative under- 
standing of the phase transition for a broad Feshbach 
resonance with a large dimensionless Yukawa coupling 
as relevant for the present experiments in ^Li and 
^''K. For zero temperature our calculations agree well 
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with Quantum Monte Carlo simulations [37| at the reso- 
nance. 

This paper is organized as follows: 

In sect. |lT]we investigate the Feshbach resonance and 
introduce the important molecule degrees of freedom in 
terms of a di-atom or molecule field (/>(a;). This al- 
lows us to cover the whole range of temperature and 
the crossover. The Bose-Einstein condensate or the su- 
perfluid order parameter corresponds to a nonvanishing 
expectation value {(j)). We establish the equivalence of 
our formulation with a purely fermionic formulation for 
which the effective interaction between the atoms con- 
tains a nonlocal piece reflecting the molecule exchange. 
Only in the broad resonance limit /i^ — > oo this interac- 
tion becomes pointlike. 

Our functional integral formulation in terms of an in- 
dependent field 4>{^) is particularly well adapted to the 
crossover from a BEG to a BCS condensate: Molecules 
and Cooper pairs are described by the same field. Of 
course, the dynamical properties depend strongly on the 
BEC or BCS regime. In particular, we compute in sect. 
mil the gradient coefficient A^f, which determines the gra- 
dient contribution to the free energy for a spatially vary- 
ing molecule field. For the BEC regime, is dominated 
by the "classical value" , corresponding to the dominance 
of the tightly bound molecules. In contrast, for the BCS 
regime the fiuctuation effects dominate A^. In this case 
(t>{x) can be associated to a collective degree of freedom 
(Cooper pair). One could omit the classical contribution 
to A^ such that (j) becomes on the microscopic level an 
auxiliary field. Indeed the presence of a molecular bound 
state is no longer crucial in the BCS regime. One may 
work in a purely fermionic setting with a local interaction 
which depends on a magnetic field B. For a broad Fesh- 
bach resonance this feature holds for the entire crossover 
region. 

We discuss the derivative expansion of the effective ac- 
tion in sect. IIIII In particular, sect. IIII Al addresses 
the issue of additive renormalization of the detuning and 
sect. IIII B I computes the wave function renormalization 
which distinguishes the "renormalized field" for the 
dressed molecules and the field for the microscopic of 
"bare" molecules. We turn in sect. IIVI to the discussion 
of the relevant parameters that describe the universal 
aspects of ultracold atoms. In terms of the dimension- 
less concentration c and Yukawa coupling the system 
becomes independent of the detailed short distance prop- 
erties. 

Sect. |V| discusses the fraction of atoms bound in 
molecules. It is crucial to distinguish between the "bare" 
or microscopic molecules and the dressed molecules 
(3S . i39j . Their numbers are related by the multiplicative 
wave function renormalization Z^. In order to achieve a 
complete symmetry between the fermionic fiuctuations of 
unbound atoms and the bosonic fluctuations of molecules 
we adapt our functional integral setting in sect. IVII In 
sect. I VIII we turn to the BEC limit. For a broad Feshbach 
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resonance (large h^) one finds very large Z^, such that 
the condensate fraction (condensed dressed molecules) 
exceeds by far the number of microscopic molecules. The 
latter becomes completely negligible for oo. Nev- 

ertheless, we find a Bogoliubov theory for bosons in the 
low temperature BEC regime for all values of /i^. Fi- 
nally, we include the molecule fluctuations (or collective 
fluctuations of di-atom states) in sect. IVIIII in the form 
of new bosonic gap equations. We draw conclusions in 
sect. US 

While the main part of this paper deals with a ho- 
mogeneous situation our formalism can be extended to 
cover the inhomogeneous situation in a trap of atoms if 
the inhomogeneity is not too large. Since the main part 
is independent of the discussion of inhomogeneities we 
display the formalism for inhomogeneous situations in 
appendix |^ 

We introduce a general formalism for a functional in- 
tegral which applies to arbitrary fermionic systems and 
is easily generalized to systems with bosons, far beyond 
the particular case of a Feshbach resonance (where di- 
atom states play a role). In addition to the (Grassmann) 
field variables ip{x) for the fermionic atoms we employ 
a bosonic field variable (t{x). It corresponds to a vary- 
ing effective chemical potential which is associated to the 
density field n(x). This procedure allows computations 
for the inhomogeneous setting of atoms in a trap beyond 
the small density approximation or beyond the Thomas- 
Fermi approximation. The bosonic field variable cr is 
introduced by partial bosonization. We formulate the 
effective action r[a] and establish the exact formal rela- 
tions between a{x), n{x), the chemical potential fi and 
the local trap potential Vi (x) . 

II. DILUTE GAS OF ULTRACOLD ATOMS 

The ultracold gas of fermionic atoms in the vicinity of 
a Feshbach resonance can be treated in the idealization 
of two stable atomic states denoted by a two component 
spinor i/j. (For the example of ^Li these states may be 
associated with the two lowest hyperfine states |1), |2).) 
The molecular state responsible for the Feshbach reso- 
nance can be treated as a bosonic particle. In our ide- 
alization it is stable for negative binding energy and can 
decay into a pair of fermionic atoms for positive binding 
energy. 

For a realistic description our formalism has to be ca- 
pable to describe the following phenomena: (i) Conden- 
sates of atom pairs may form at low temperature, sim- 
ilar to the ECS description of superconductivity, (ii) 
Molecules of two atoms can be exchanged between the 
single atoms, thus contributing to the interaction. Also 
these molecules may form a Bose-Einstein condensate at 
low temperature. In our formalism both effects find a 
unified description as will be discussed in detail in this 
paper. 

We work with a microscopic action which explicitly 



includes a bosonic field (f> with atom number two (25j . 

Sb ^ J dx[^^idr-^A-a)i^ 

-^h{r^'^ei,-^i^hr)}. (1) 

Here we employ the Matsubara formalism where the Eu- 
clidean time r is wrapped on a torus with circumference 
f3 — 1/T with conventions 

a;=(T,x), Jdx^J dr J d^x. (2) 

(Our units are h = c = = 1-) We will shortly see how 
this microscopic model relates to a purely fermionic set- 
ting by the means of a partial bosonization or Hubbard- 
Stratonovich transformation. 

The complex two-component spinors tpi^) anticom- 
muting Grassmann variables. We assume an equal mix- 
ture of the two atomic states. In this case the chemi- 
cal potential associated to the difference in the number 
of atoms in the "up" and "down" states precisely can- 
cels the energy difference between the two states such 
that both can be omitted. (For unequal mixtures the ac- 
tion contains an additional term cx ip'^r^il).) The bosonic 
molecules are described by a complex bosonic field 0. 
The propagator for these "bare molecules" is obtained 
from simple symmetry considerations, i.e. we assume 
a mass 2M, leading to a nonrelativistic kinetic energy 

/AM . The quadratic term ~ (jfcj) involves the "bare" 
binding energy or detuning {va) which typically depends 
on the magnetic field. In order to make contact to phys- 
ical observables, has to be additively renormalized, 
which is implemented in sect. IIII Al 

We use here two different chemical potentials a and /i 
for the fermionic atoms and bare molecules. This is useful 
if we want to obtain separately the densities of fermionic 
atoms or bare molecules by differentiation of the free en- 
ergy with respect to a or /i. Since only the total number 
of atoms is conserved, one has to set a = fi &t the end 
of the computations. The distinction between a and /i is 
appropriate if one wants to understand explicitly the role 
of the microscopic (or bare) molecules. In sect. IVII we will 
drop this distinction in favor of a more unified approach. 
There we will set ii — a from the outset such that a will 
be conjugate to the total density of atoms irrespective 
of a microscopic distinction between unbound atoms and 
molecules. In the main part of this paper we will treat 
a and /i as constant classical source terms. However, we 
stress that this source term can be straightforwardly pro- 
moted to a fiuctuating field, a — > a{x). This issue, and 
its use for the description of inhomogeneities beyond the 
usual Local Density Approximation, will be discussed in 
the appendix [X\ 

The Yukawa or Feshbach coupling describes the 
coupling between the single atoms and molecules, eap = 
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— e/3a, £12 = 1- For ^ the molecular states decouple. 
However, in an appropriately performed "narrow reso- 
nance limit" /i0 — > which keeps the scattering length 
fixed, an exact solution of the many-body problem be- 
comes feasible above the critical temperature. A detailed 
analysis of this limit is given in [3] ■ Broad Feshbach res- 
onances correspond to large h^, and we will see that the 
limit oo describes a purely fermionic theory with 

pointlike interactions where microscopic molecules can be 
neglected. 

The thermodynamic equilibrium situation is described 
by the partition function. The basic ingredient for our 
formalism is the representation of this object in terms of 
a functional integral with weight factor e^"^^, with Sb 
the Euclidean action ^ 

Zhj] = j V^V^e^^[-SB[iPA] (3) 

-|- Jdx ri{x)'iip'^ (x) + rj^ip{x) 

+j{x)$*{x)+j*{x)${x)y 

This formulates the full quantum theory in terms of the 
sources ry and j for the fermion and the boson fields. 
All one particle irreducible n-point functions, including 
the order parameter and the correlation functions, can 
be directly extracted from the effective action F, which 
obtains by a Legendre transform, 

T[iP,4>] = -lnZ + Jdx T]{x)tp^ (x) + T]^x)ip{x) (4) 

+jix)^*{x) +j*{x)^ix). 

The effective action is a functional of the "classical" fields 
(or field expectation values) tp = {'ip) ,4> = {(p) in the 
presence of sources. They are defined as 

0(x) = (0)(x) = |^ (5) 

and analogous for the fermion fields. Of course, due to 
Pauli's principle, the fermion field cannot acquire a non- 
vanishing expectation value for rj = rj^ — such that the 
physical value is simply tp = 0. It is often convenient to 
write F as an implicit functional integral over fluctuations 
Sep, 5ip around "background fields" (p, ip 

T[i^, ^] = - In j VSi:V6<pexp {~S[ij + Sip, ^ + 5(P] + 

J if S(P + r]^^ + h.c), (6) 

with j*{x) = ST/S(p{x). This form is particularly useful 
for the construction of the equation of state, i.e. the 
explicit expression for the total atom number density. 

The total number density of atoms n includes those 
from unbound or "open channel" atoms and the ones aris- 
ing from the bare molecules or "closed channel" atoms. 
It obeys 

n{x) = npix) + neix) = {ip\x)ijj{x)) + 2{$*{x)4>{x)).{7) 



Indeed, the action is invariant under t/(l) phase trans- 
formations of the fermions and bosons, 

^ ^ e"'V', ^ ^ e^'U (8) 

and the corresponding Noether charge is the total atom 
number N = j dPxn{x). We emphasize that eq. ([7|) is 
no "ad hoc" assumption - it is an exact expression for 
the particle number and directly follows from the mi- 
croscopic formulation. More technically speaking, (ip^ip) 
and 2{(p*(p) represent the full two-point correlation func- 
tions of the "bare fields" which appear in the microscopic 
action ([T]) and are quantized by means of the functional 
integral. In a homogeneous situation, the conserved par- 
ticle number can be replaced by a fixed constant particle 
density, n — N/V. 

The bosonic part n b counts the total number of atoms 
contained in the microscopic or "bare" molecules. This 
number receives a contribution from free molecules and 
from the condensate as discussed in more detail in sect. 
IVl In the language often used for a Feshbach resonance, 
ub measures the "closed channel" microscopic atoms. 
Using the formalism of the present paper we have com- 
puted fiB as a function of the magnetic field B in p^ . 
We find very good agreement with observation [§| over 
several orders of magnitude in hb- 

Since the action ([1]) contains only terms quadratic or 
linear in (p it is straightforward to express the expectation 
value (pQ in terms of a local fermion-bilinear. It obeys (for 
vanishing source for (p) 

(i^A - 2/i - ^ + dr)^0 - ^(V^^eV^). (9) 

In particular, for constant (pQ we find 

^o^-^i^ei^). (10) 
zi^A — 4/x 

This demonstrates directly that our formalism makes no 
difference between a "condensate of molecules" (p and 
a "condensate of atom pairs" {iP^eip) - they are simply 
related by a multiplicative constant. 

We finally show the equivalence of our formalism with 
a model containing only fermionic atoms and no micro- 
scopic molecules. The interaction in this fermionic de- 
scription is, in general, not local. It becomes local, how- 
ever, in the limit of a broad Feshbach resonance for /i^ — > 
oo. For this purpose we use again the quadratic form of 
the bosonic part of the microscopic action ([T]), which al- 
lows us to integrate out the (p field. Expressed only in 
terms of fermions our model contains now a momentum 
dependent four- fermion interaction ((54 = Qi-\- Q2 — Q3', 
Q = (ujn, q) with discrete bosonic Matsubara frequencies 
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ll}„ = 2'KnT at finite temperature) 

s^nt ^ -\ I {i^\-Qi)i>{Q2)){i^\QM{~Q-i)) 



I Df, - 2/i + 27ri(ni - n^)T + {qi - qi)^/4:M J ' 

(11) 

We emphasize that there is no diff'erence between the 
Yukawa type model described by the action ([T]) and a 
purely fermionic model with interaction (jlip . All physi- 
cal observables can be computed in either one or the other 
of the two formulations. However, eq. (jlip reveals that 
our model describes a setting beyond pointlikc interac- 
tions via the classical frequency and momentum depen- 
dence of the four-fermion interaction. The momentum 
structure of (fTT|) is compatible with interactions in the 
tpip - channel. The action ([T]) hence models a nonlocal 
coupling between the fermionic constituents. Reversing 
the logic, eq. (H]) could also be obtained by starting from 
eq. (jlip and performing a Hubbard-Stratonovich trans- 
form or partial bosonization [13, |4l|. Finally, we note 
that we could also choose classical "gradient coeSicients" 
A^"'-* (see below) different from 1 / (4M) in order to model 
an experimentally determined effective range. 

In the pointlike limit the momentum dependence and 
the dependence on ^ can be neglected and the coupling 
term in eq. (fTTj) is replaced by the "local interaction 
approximation" 

Aa = (12) 

This limit obtains formally for h'^ —> oo, D/^ —> oo, while 
keeping A fixed. It is relevant for broad Feshbach reso- 
nances, as discussed in detail in (19] . 



III. DERIVATIVE EXPANSION FOR THE 
EFFECTIVE ACTION 

The condensation of atoms pairs or molecules is sig- 
nalled by a non- vanishing expectation value {(f)) = 4>o. 
The associated symmetry breaking of the global contin- 
uous symmetry of phase rotations of ip and (p (related 
to the conservation of the number of atoms) induces a 
massless Goldstone boson. This is the origin of super- 
fiuidity. In this section, we show how to describe this 
phenomenon in the effective action formalism. For this 
conceptual issue, it is sufficient to work in the mean field 
approximation or a simple extension thereof (extended 
MFT). The more sophisticated approximation schemes 
beyond mean field are presented in sects. IVI[ IVIIII 

In this work we will treat the effective action in a 
derivative expansion, i.e. we write 



and compute the "wave function renormalization" Z^, 
the "gradient coefficient" and the effective potential 
U{(f>). In the present paper we do not compute the correc- 
tions to the part of the effective action involving fermions 
- for this part we have simply taken the classical action. 
(See ll^ for the renormalization of the Yukawa coupling 
in presence of a "background" four-fermion interaction.) 
We therefore omit the fermionic part of the effective ac- 
tion from now on. For the concrete calculation we work 
in momentum space. We emphasize, however, that the 
above expression can be used for the investigation of weak 
inhomogeneities as encountered in atom traps. The fluc- 
tuation problem, i.e. the computation of Z^,A(f,,U in 
the above truncation, can then be solved in momentum 
space, while the effects of weak inhomogeneities can be 
investigated by solving the classical field equations de- 
rived from (fT3)) . This reaches substantially beyond the 
usual local density approximation, which ignores the ki- 
netic terms in eq. (|13p . We discuss the implementation 
of an external trapping potential in our functional inte- 
gral formalism in app. El Here, however, we focus on 
the homogeneous situation. 

In its most general form, the effective action depends 
on the parameters T, a and ^ and on the classical field (j). 
The dependence on the chemical potentials a and ^ can 
already be inferred from the partition function - they are 
only spectators w.r.t. the Legendre transform. Follow- 
ing the thermodynamic construction, the total particle 
density in a homogeneous setting is obtained as 



dU_ 
da 



(14) 



This prescription precisely reproduces eq. Here U 

has to be taken at its minimum. 

In the absence of sources the field equation for ^(x) 
reads 



ST 
54>{x) 



= 0. 



(15) 



For a homogeneous situation the stable solution corre- 
sponds to a minimum of the effective potential 



with 



op 



(16) 



(17) 



Here we have introduced the U{1) invariant p ~ <j)*<f> - 
the effective potential only depends on this combination. 
This simple equation can be used to classify the thermo- 
dynamic phases of the system. 



— dx 



{i7(0) + Z^^*dr4> + A^v^*y^ + ...},(13) 



Symmetric phase(SYM) 
Symmetry broken phase(SSB) 
Phase transition(PT) 



Po = 0, 
Po > 0, 
Po = 0, 



U'iO) > 0, 
U'{po) = 0, 
[/'(O) = 

(18) 
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where po = (I)q4>o denotes the minimum of U{p). For high 
temperatures, the minimum of the effective potential oc- 
curs at = and we deal with a normal gas phase. For 
low enough T, on the other hand, we expect the mini- 
mum of U{a, 4>) to occur for 0o 7^ 0- The spontaneous 
breaking of the U{1) symmetry is signalled by a nonzero 
field expectation value and indicates the condensation 
phenomenon. This has an important aspect of universal- 
ity: one and the same criterion can be used for the whole 
parameter space, both for the BCS-type condensation of 
Cooper pairs and for the BEC of microscopic molecules. 

In the remainder of this section, we will evaluate the 
effective action in the mean field approximation (MFT). 
This scheme is defined by only considering the effects gen- 
erated by fcrmion fluctuations. We will include bosonic 
fluctuations in later chapters. Beyond MFT, we first 
include the contribution to the density from dressed 
molecules. This is connected to the effective bosonic 
two-point function (connected part of the bosonic par- 
ticle density). This effect is included in different current 
approaches to the crossover problem in the limit of broad 
Feshbach resonances [l^, [23, HI] . We will refer to it as 
extended Mean Field Theory. Furthermore, we include in 
sect. I Villi the modifications of the effective potential due 
to bosonic fluctuations, using suitable Schwinger-Dyson 
equations. 

In a realistic physical situation the validity of our 
model is restricted to momenta smaller than some mi- 
crophysical "ultraviolet cutoff" A. In turn, A is typi- 
cally given by the range of the van der Waals interac- 
tions. One may use A « as/lOO with ag the Bohr ra- 
dius. For practical computations it is often convenient to 
consider the limit A —^ 00 such that no explicit informa- 
tion about the "cutoff physics" is needed. This requires 
to express the couplings of the theory in terms of suit- 
ably defined "renormalized couplings" that stay finite for 
A — > 00. In the next subsection we will discuss the ad- 
ditive renormalization of the detuning. The functions 
Zcf,,A^, in contrast, are UV finite. We will also subtract 
field-independent pieces linear in a and n that obtain in 
a naive MFT computation. This additive "density renor- 
malization" will be traced back to the relation between 
the functional integral and the operator formalism. 



A. Effective potential and additive renormalization 

In the mean field approximation, the effective potential 
reads, after carrying out the Matsubara summation and 
omitting an irrelevant infinite constant, 

Ua{<j,^) = (j>A-2/i)0*0 + AC/|^\ (19) 
^Uf^ = -2T / — ^lncosh7^ 



where 




Here we have added an index A in order to remind that 
this form still depends on an ultraviolet cutoff A. We reg- 
ularize the effective potential by limiting the integration 
over spacelike momenta by an upper bound A, (/■^ < A^. 

Let us now discuss the additive renormalization needed 
to properly describe the physics encoded in this object. It 
is needed in two instances: The first one concerns a zero- 
point shift of the two-point function and is related to 
the quantization via the functional integral. Its removal 
does not involve a physical scale and can thus be seen 
as a normalization of a certain observable, the particle 
number. The second one is related to a true ultraviolet 
divergence which needs to be cured by an appropriate 
counterterm. For the ultraviolet renormalization we can 
restrict to the simpler situation where there is no spon- 
taneous symmetry breaking, since this effect occurs only 
in the low energy sector and cannot affect the ultraviolet 
physics. For the physical value of the field expectation 
value, this implies <^ = and 7^ = 7. 

The fermionic part of the particle number is naively ob- 
tained from the thermodynamic relation ua = —dUt^/da 
in a homogeneous setting. This yields the explicit ex- 
pression 

nF.A = - / 7-— ^tanh7 (21) 

Ja i^T^r 

and we observe that this number may get negative for 
large negative a/T. In order to clarify the precise rela- 
tion between h\ and the particle density, we first consider 
the simpler situation of a single fermionic degree of free- 
dom. In this case the expectation value can be 
related to the expectation values of products of the usual 
annihilation and creation operators a, a^, which obey the 
anticommutation relation a^a + aa^ ~ 1, 

{ijj'^iil)) = — (V^^Vi — = — (a^a — aa^) 

= (a^a) - 1/2 = n - 1/2. (22) 

Here the second equality holds since this combination of 
operators is covariant with respect to permutations of 
the ordering. For a lattice model (as, for example, the 
Hubbard model) with / degrees of freedom per site the 
fermion number per site therefore reads n = n\ + ^. 
For electrons in a solid (/ = 2) one can associate n\ 
with the difference of the electron density from half filling 
(where n = 1), i.e. the average number of electrons minus 
holes per site as compared to the half-filling density. For 
relativistic charged fermions '^''''ijj measures the difference 
between particle and antiparticle density and the additive 
constant drops out. For nonrelativistic atoms, however. 
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the relation between the atom density n and n\ becomes 

3 

UfA^) = {■4>Hx)iPix)) (23) 

= ^(/E [^1(^)V^j(2/) - My)^}i^)] - y)) 

V 
y 

= (a^(a;)a(x)) - -(5(0) = ^^(a;) - - ^ ^^-p- 
= np{x) — n. 

The volume factor in momentum space, 5(0), diverges 
in the limit of infinite momentum cutoff. The physical 
fermionic particle density 71^(2;) and the relative particle 
density np^x) = {tp"^ {x)ip{x)) are therefore related by 
an additive shift that depends on the momentum cutoff. 

In consequence, one finds now a manifestly positive 
total fermionic particle number 

Nf = j d?x{nF,K + n)^V j ( exp(27) + l) 

(24) 

The momentum integral is now ultraviolet finite for A — > 
00. It becomes exponentially insensitive to the ultravio- 
let cutoff A, such that we dropped the index. We recover 
the familiar Fermi distribution. An analogous argument 
holds for the connected part of the bosonic two-point 
function which will be implemented below. Formally, 
this additive renormalization appears in the form of field 
independent terms linear in a and fi in the classical po- 
tential. 

Let us now proceed to the second instance where UV 
renormalization is needed. The microscopic action p]) 
depends explicitly on two parameters v\, h^. A third 
parameter is introduced implicitly by the ultraviolet cut- 
off A for the momentum integration in the fluctuation 
effects. (Besides this, the results will depend on the 
thermodynamic variables T and cr, //.) Contact to ex- 
periment is established by relating the microscopic pa- 
rameters to observables of the concrete atomic system. 
Here we choose these parameters to be the magnetic field 
dependent physical detuning v{B) and the Feshbach cou- 
pling The parameters (P, h^) can, however, be re- 
placed by an equivalent set {a~^,h^) (a the scattering 



^ Note that the constant shift h diverges if the ultraviolet cutoff 
for the momentum integration (g^ < A^) goes to infinity, h = 
A3/(67r2). 

* For vanishing background or open channel coupling, is a free 
parameter. Else, a further UV renormalization of the Feshbach 
coupling is necessary O . 



length) in a second step as pointed out below. Once 
the parameters v{B) and /i^ are fixed by the properties 
of the molecules or atom scattering in empty space we 
can proceed to compute the properties of the atom gas 
at nonzero temperature and density without further free 
parameters. 

In the vicinity of the Feshbach resonance at B = Bq 
we may approximate v\{B) by a linear behavior (linear 
Zeeman effect) 



duA 
dB 



MB- 



(25) 



Here pB — M-i- + M- ~ MAf reflects the difference between 
the sum of the magnetic moments of the two atomic 
states (p^+ii-) and the molecule magnetic moment ^m- 
We relate the physical detuning D to Pa by an additive, 
S-independent shift 



/i^MA 



di^A 
dB 



dv 
'dB' 



(26) 



This is motivated by a consideration of the fermionic con- 
tribution to the "boson mass term" , 



= C/'(po = 0). 



(27) 



Indeed, the fluctuation contribution diverges for K —t 00 



ma 



2 J {2TTf^2M 



2 



d3 



q 



{2n) 



^2M 



-tanh^!/^^ 
2T 

-tanh.'^^/^^^-'^ 



2T 



2M 



(28) 



In the second equation the linear dependence of the fluc- 
tuation correction on the cutoff A is absorbed into the 
definition of the physical detuning D. The remaining in- 
tegral in the second line of ([28|) is very insensitive with 
respect to the precise value of A, and we can formally 
send A to infinity. 

More generally, we choose vk — v such that the zero of 
V coincides with the Feshbach resonance at Bq 



V = [Ib{B - Bq). 



(29) 



In vacuum (n — 0,T — 0) the Feshbach resonance cor- 
responds to a vanishing binding energy eM = 0. This is 
realized I5| for cr = 0, fn^ = 0. There is also no conden- 



sate in the vacuum, i.e. 4>o — 0. We note that there are 
no boson fluctuations contributing to the renormalization 
of the mass term in the physical vacuum n = T = as 
can be seen from a diagrammatic argument [4^ . 

The scattering length a can be defined in terms of the 
scattering amplitude at zero momentum and zero energy 
At the present stage we may consider a "resonant 
scattering length" related to D by 



M 
Airar- 



V 



"hi 



MA 
2^' 



(30) 
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It accounts for the contribution of the molecule exchange 
to the atom scattering. For B ^ Bq and low enough 
momenta and energies (low enough temperature) the 
molecule exchange can be described by an effective point- 
like interaction. We therefore also define the renormal- 
ized resonant four-fcrmion vertex Xr, 



Comparison with eq. 



1 



4:TTaR 
M ' 

ni) yields 

Pa ma 



1 



MA 
2^' 



(31) 



(32) 



This demonstrates that the renormalization of D ([26)1 cor- 
responds to the renormalization of the effective atom in- 
teraction strength A. 

It is instructive to consider in the mean field framework 
the explicit equation which determines the order param- 
eter 00 ill the superfluid phase. A nonzero 0o obeys 



2^ 



AT J 



(27r)3 



(7^^anh7^-4MT/g2). (33) 



In an alternative purely fermionic formulation we may 
compute 00 in the local interaction approximation (jl2[) 
by solving the Schwinger-Dyson equation [i^, [4^. Eq. 
([33| would correspond precisely to the BCS gap equation 
(lowest order Schwinger-Dyson equation) for the purely 
fermionic formulation with local interaction (45l . |46| | , pro- 
vided we choose 



1 i^A — 2/i 



1 v-2n 

<t> 



(34) 



This suggests the definition of a /i or density depen- 
dent effective coupling A(^) (cf. eq. pTjl ) for the 
atom-molecule model. In the many-body context and in 
thermodynamic equilibrium, where fermions and bosons 
share a common chemical potential tr = /i, the latter is 
determined by the density. In the physical vacuum, ob- 
tained by sending the density to zero, describes the 
binding energy of a molecule, cm = 2/i = 2(T and does 
not vanish on the BEC side of the resonance [IB]- We 
emphasize that the resonant scattering length Xr in eq. 
(Plj) describes the scattering of fermions throughout the 
crossover and is directly related to the observed scat- 
tering length for two atom scattering. On the other 
hand, A(/x) is a universal combination characteristic for 
the ground state of the system 42] . For broad Feshbach 
resonances (/i^ oo) the two quantities coincide. 

Expressing the effective potential U (|19p in terms of v 
the momentum integral becomes ultraviolet finite. 



-2T 



(27r)3 



ln(e'^*"^ -f e"'^'*-'') 



2Tg2 



(35) 



The remaining cutoff dependence is 0{A~^), the precise 
value of A therefore being unimportant. For definiteness, 
the cutoff A can be taken A = {Sn'^ngsy^^ with Ugs the 
density in the liquid or solid ground state at T = 0. This 
choice will be motivated in app. [B] 

At this point we have reached a simple but nevertheless 
remarkable result. When expressed in terms of measur- 
able quantities, namely scattering length (or Xr) and 
the open channel atom density np the (p dependence of 
the effective potential becomes very insensitive to the mi- 
croscopic physics, i.e. the value of the cutoff A. Without 
much loss of accuracy we can take the limit A — > cxd for 
the computation of U. The effective chemical potential 
a depends on np and (p via dU{a,(l))/da — — fip- 



B. Wave function renormalization 

The wave function renormalization is another im- 
portant ingredient for the description of the crossover 
physics in the Yukawa model. This can be seen from the 
fact that rescaling all couplings in the effective action 
with the appropriate power of .Z^, we end up with an 
effective bosonic Bogoliubov theory in the BEC regime, 
as will be discussed in detail in sect. IVIII Here we focus 
on the explicit computation of the wave function renor- 
malization. 

As can be read off from eq. (fT5|) . the wave func- 
tion renormalization is related to the inverse molecule 
propagator 7^0. In app. [C]we compute 7^^ in the mean 
field approximation. 



V^{Q) = 27rinT+^ 



(36) 



PF{Q')-'y 



Pf{-Q' + Q)-(t 



[Pf[Q' - Q) - a][PF{-Q' + Q)-a] + hl4>*4> 

where the kinetic part of the inverse fermion propagator 
reads 



q 

2M 



(37) 



with Q — (ujB,q), and the frequency variable lob rep- 
resents the discrete fermionic Matsubara frequencies at 
finite temperature, ujb — (2n + 1)ttT. 

Here we interpret the wave function renormalization 
^ as a renormalization of the term in the effective ac- 
tion with a timelike derivative. We may then evaluate the 
propagator correction (|36p for analytically continued fre- 
quencies ll>b ^ + it^ and set ujb = 0. Now lS.P^{u!, q) 
becomes a continuous function of lo. Defining 



Z^ = 1 



dAP^ 
du! 



(38) 
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one finds 



16f 2 



(^77^ 3 [ tanh 7^ - 7^ cosh ^ 



70J 
(39) 

Here we have rescaled the integration variable and the 
Yukawa coupling as g = g/^F, = AM'^h'^/kp. In the 
symmetric phase the simplification 7^ = 7 applies. We 
note that is closely related to the spectral function for 
the molecules. If we only consider fermionic diagrams, 
we can give an equivalent definition of the wave function 
renormalization using the mean field effective potential, 

1 d^U 

Zl = \-- - - . (40) 
^ 2dad(t>*d(l) ^ ' 

The property — Z^ holds since the integral in eq. 
pop depends on the combination uj — 2(t only. 

C. Gradient coefficient 

In order to compute the gradient coefficient, we pro- 
ceed in complete analogy to the wave function renormal- 
ization: For the spacelike momenta, we define 



Ml) 



^40, g)- ^^(0,0) 



lim A^{q). 



More explicit formulae are given in app. 
(I013I014|) . 



(41) 



[Cl in eqs. 



IV. RELEVANT PARAMETERS, MOMENTUM 
AND ENERGY SCALES 

A. Concentration 

For a mean field computation with nonzero density the 
effective renormalized four-fermion coupling is related by 
an equation similar to eq. (|12p to the resonant scattering 
length in vacuum. We define ^ 



1 



M 



A(cr) 47ra_R 



2a 



(42) 



Therefore the effective scattering length in an atom gas 
differs from the (vacuum) scattering length which is mea- 
sured by the scattering of individual atoms. This den- 
sity effect is refiected by the dependence on the effective 
chemical potential and we define 

1 _ 1 Sttct 
a gr 



hlM' 



(43) 
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FIG. 1: Crossover phase diagram ^^(c"^). For large we 
compare a calculation with a gap equation modified by boson 
fluctuations (dashed) to the result obtained with the standard 
BCS gap equation (solid). The universal narrow resonance 
limit /i0 ^ is indicated by the dashed-dotted line. We fur- 
ther plot the result of the standard BCS approach (BCS gap 
equation, only fermionic contributions in the density equa- 
tion; long-dashed rising line) and the result for noninteract- 
ing bosons (long dashed horizontal line). The solid line at 
resonance Tc{c~^ = 0) = 0.292 coincides with the result ob- 
tained in which omits boson fluctuations. The dashed line 
(T'c(c~^ = 0) = 0.259 is in good agreement with the value ob- 
tained in [l^l (Tc = 0.266), who work in the broad resonance 
limit from the outset. 



On the BCS side a is positive and is negative - the 
size of |a| is therefore larger than |afl|. A similar enhance- 
ment occurs in the BEC regime where a will turn out to 
be negative and positive. Roughly speaking, on the 
BCS side the presence of a nonvanishing atom density 
favors the formation of (virtual) molecules by reducing 
the "cost" of forming molecules with positive v (or pos- 
itive binding energy) to v — 2a. In the BEC regime the 
presence of molecules reduces the absolute size of the ef- 
fective binding energy. It should, however, be pointed 
out that the effect of a density dependence is small in 
case of broad resonances — s- 00, as visible from eq. 
dMl). 

The atom density n defines a characteristic momentum 
scale by the Fermi momentum kp, i.e. 



l3 

3^' 



(44) 



The inverse of the Fermi momentum defines the most 
important characteristic length scale in our problem. 
Roughly speaking, it corresponds to the average distance 
d between two unbound atoms or molecules. We empha- 
size that our definition of kp involves the total density 
n and therefore includes unbound atoms, molecules and 
condensed atom pairs. In terms of kp we can form a 
characteristic dimensionless "concentration" 



For a = fj, this combination appears in the term quadratic in 
in eq. I I35I I. 



afci 



(45) 



10 



The concentration is a measure for the ratio between 
the in-medium scattering length a and average distance, 
c (x a/d. As mentioned in the introduction, the concen- 
tration is the crucial parameter for the description of the 
crossover between the BEC and BCS regimes. For a small 
concentration |c| the gas is dilute in the sense that scat- 
tering can be treated as a perturbation. In this range 
mean field theory is expected to work reasonably well. 
On the other hand, for large |c| the scattering length ex- 
ceeds the average distance between two atoms and fluc- 
tuation effects beyond mean field might play a crucial 
role. 

An alternative definition of the concentration could use 
the measured (vacuum) scattering length, c = a{B)kF- 
This definition has the advantage of a simple relation to 
the magnetic field. The choice between the two defini- 
tions is a matter of convenience. We adopt here the def- 
inition (j45p since this refiects universality in an optimal 
way. For broad Fcshbach resonances the two definitions 
coincide since a — a{B). In presence of an additional 
"background scattering length", a{B) — ai,g + aji{B) we 
will include a^g in the definition of c. 

The concentration c is the most important parame- 
ter for the description of the crossover (besides T and 
n). The inverse concentration corresponds to a "rel- 
evant parameter" which vanishes at the location of the 
resonance. Once described in terms of c^^ the ultracold 
fermionic gases show a large universality. We demon- 
strate this universality on the crossover phase diagram in 
fig. [T]which plots the critical temperature Tc = Tc/ep for 
the transition to superfiuidity as a function of c^^. The 
narrow resonance limit is exact, while our results 

for broad resonances {h^ — > oo) still have substantial un- 
certainties, as demonstrated by two approximations that 
will be explained later. For large the actual value of 
becomes irrelevant and all curves coincide with the 
broad resonance limit. Intermediate /i^ interpolate be- 
tween the broad and narrow resonance limits. 

We have argued in sect. |TT]that in the limit of a point- 
like approximation for the effective fermionic interaction 
all results should only depend on the effective scatter- 
ing length. This only involves the ratio i^/h^ such that 
for fixed c the separate value of should not matter. 
On the BCS side for small |c| we therefore expect a uni- 
versal behavior independent of the value of the Yukawa 
coupling h^. This is clearly seen in fig. [1] where the crit- 
ical line approaches the BCS result independently of h^. 
Furthermore, all results become independent of the value 
of in the broad resonance limit {h^ —^oo). The con- 
centration remains then as the only parameter (besides 
T and n) for the description of the crossover. These new 
universal aspects, adding to those that will be presented 
in IIV C[ are discussed in [l^. In particular, the broad 
resonance limit /i^ — > c» corresponds to a pointlikc mi- 
croscopic interaction. 

For the broad resonance limit the first approximation 
(solid line) corresponds to extended MFT and neglects 
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FIG. 2; Gradient coefficient = 2MA^ and wave function 
renormalization in dependence on . We divide by h?^ 
in order to get numbers 0(1) and use — 3.72 ■ 10^ as 
appropriate for ^Li p^ . The ratio — A^/Z^ is displayed 
in fig. H 

the modifications of the effective potential induced by 
the molecule fluctuations. The second approximation 
(dashed line) includes these fluctuation effects via the 
solution of a gap equation for the molecule propaga- 
tor (sect. IVIIip . The fast approach to the BEC value 
for 0+ does not reflect the expected behavior 

Tc = T^^'-^ + Kc with a dimensionless constant k. This 
shortcoming of our treatment is remedied by a functional 
renormalization group study [49| . The solid line at res- 
onance Tc{c~^ = 0) = 0.292 coincides with the result 
obtained in [13] which omits boson fluctuations. The 
dashed line (Tc(c~^ = 0) = 0.259 is in good agreement 
with the value obtained in ^] {Tc — 0.266), who work 
in the broad resonance limit from the outset. 

In fig. [2]we show the wave function renormalization 
and the dimensionless gradient coefficient = 2MA^ 
as a function of c^^ for T = 0. We use the large value for 
the Feshbach couphng in ^Li, = 3.72-105 for kp = leV 

(cf. eq. (|46|) ). We note that for large h^, as appropriate 
for the Feshbach resonances in ^Li or ^"K, the wave func- 
tion renormalization is large, the ratio Z^/h?^ being an 

0(1) quantity. We observe a strong increase of when 
evolving to the BCS side of the resonance. This strongly 
suppresses the propagation of the effective bosonic de- 
grees of freedom, leading to a situation where they are 
completely irrelevant. This is an aspect of universality, 
where the system "looses memory" of the bosonic degrees 
of freedom and a purely fermionic BCS-type description 
becomes appropriate. These curves are essentially iden- 
tical in the limit /i^ ^ cx3. 

Indeed, large values of h?^ exhibit an "enhanced uni- 
versality" [13]. In this case, all "microscopic quantities" 
depend only on the concentration c. In the limit hff, oo 
at fixed scattering length there is a loss of memory con- 
cerning the details of the bosonic sector in the micro- 
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scopic action ([T]). These aspects are worked out in more 
detail in [TH . In a renormalization group treatment this 
universahty will be reflected in strongly attractive partial 
infrared fixed points, similar to the quark meson model in 
strong interactions [50J. This universality property will 
be valid beyond mean field theory. 



approximation the effective chemical potential a obeys, 
using the dimcnsionless shorthands (|47p . 



dqq { — tanh70-l} = -- . (51) 

70 3 n 



B. Dimensionless parameters 



The characteristic energy scales for T, a and the gap 
A = h^(j) a^re set by the Fermi energy ep — kp/{2M). It 
is appropriate to define dimensionless quantities 



T 



2MT/kl = T/ep, = 2Ma/kj,, 
q/kp, A/ep ^2Mh^4>/kl, f = A*A, 

(46) 



= 2Mkp^''^h^, Aa,^2MAa 



Once all quantities are expressed in this way in units of 
the Fermi momentum or the Fermi energy the atom mass 
M will no longer be present in our problem. Indeed, in 
dimensionless units one has for the expressions (|20p 
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1 (~2 

2T^ 
1 



(47) 



70 - 



.X 1/2 

r 



such that the atom mass M drops out in the computation 
(|19I35P of the appropriately rescaled effective potential 
U . All dimensionless quantities in 7^ are typically of the 
order one. For practical computations we may therefore 
choose units kp = leV. 
The rescaled potential 



u = kp''-U ^ 2MkpW 



(48) 



is composed of a classical contribution ^ and a contribu- 
tion from the fcrmion fluctuations {u^^"^) 



(49) 



1 



fj — 



Sttc Sttc h'^Mkp ' 

Then the equation determining a{r,T) becomes 

du 1 fip 

da Stt^ n 



(50) 



and is indeed independent of M . It depends, however, 
on the ratio np/n since we have defined in eq. (j44p kp 
as a function of n rather than rip. In the mean field 



In the formulation of sect. IVII c will bo replaced by c since u and 
will be identified. 



In particular, for f = this determines (i as a function 
of T and npjn. 

For given a, T and c one can compute the order pa- 
rameter f in the low temperature phase according to 



1 



du ^ du\ 
df ' df Sttc 



(52) 



The critical temperature for the phase transition corre- 
sponds to the value where f vanishes. This part of the 
mean field computation is independent of the Yukawa 
coupling In sect. |V] we will determine fipjn as a 
function of (t, T and c such that eq. (jSOp can be used 
to determine ct as a function of c and T. As an alter- 
native, we will modify in sect. IVII the definition of V 
such that a modified equation du/da = — I/Stt^ becomes 
independent of npjn. Again, a can now be fixed as a 
function of c and T. We will see that the relation (t(c, T) 
depends on the choice of the Yukawa coupling h^. The 
results will therefore depend on the additional dimension- 
less parameter . Away from the narrow and broad 
resonance limits, the model is therefore characterized by 
two dimensionless quantities, c and h^. 



C. Universality 

All observables can be expressed in terms of c, T and 
h^. For example, using the definition 



2M _ 



one finds the relation 



v-2a 



Sttc 



(53) 



(54) 



In particular, the phase diagram Tc{c) depends only on 
as shown in fig. [T] for the case of narrow and broad 
resonances. For T — Tc find that the relation a{c) 
depends only mildly on the value of ft^ as shown in fig. 
[3] where we compare narrow and broad resonance limits 
again. Furthermore, for small |c| the universal curves 
Tc(c) and (t(c, T) become independent of h^: Both in the 
BEG and BCS regime, an enhanced universality sets in, 
making h^j, irrelevant for all its possible values! These 
issues are investigated more systematically in [isj . 

In summary we have now an effective low energy for- 
mulation where neither M nor A enter anymore. Every- 
thing is expressed in terms of kp and three dimensionless 
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FIG. 3: Crossover at the critical temperature in the broad 
resonance limit: Effective dimensionless chemical potential a 
at the critical temperature as a function of the inverse con- 
centration . We compare the results for two versions of 
the gap equation as in fig. [1] (solid and dashed line) . Addi- 
tionally, the result for the narrow resonance limit is indicated 
(dashed-dotted) . 



parameters, namely c, and T. This scaling property 
Is an important aspect of universality In H we discuss 
the relation of the parameters c and with the physi- 
cal observables for an experimental setting, namely the 
magnetic field B and the binding energy in vacuum. In 
the present paper we treat as a free parameter. The 
values of for ^Li and ''"K turn out to be large such 
that the broad resonance limit applies to these systems. 

We display in tab. [J the values of the dimensionless 
scale rations M/kp, A/kp and ep/kp. The first two give 
an idea how well the detailed microphysics decouples for 
realistic ultracold fermionic gases. We use a density n — 
4.4 • IQi^cm-^, kp = leV=(3.7290 • lO^as)"^ 



TABLE I: Typical values for the dimensionless scale ratios for 
•^Li and ''"K {kp = leV=(3.7290 ■ lO^as)"^). 





M/kp 


A/kp 


ep/kp 


•^Li 


5.65 • 10^ 


1.6 ■ 10^ 


8.9- 10"" 


«K 


40.0 • 10^ 


1.2 • 10^ 


1.2 • 10~" 
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FIG. 4: Momentum dependence of the gradient coefficient in 
the broad resonance regime > 1. We plot = 2MA^/Z4, 
for T = Tc in the different regimes, BCS (c^^ = -1.5, = 
0.057, topmost), crossover (c"^ =0,1^ = 0.29) and BEG 
c'^ = 1.5, fc = 0.218). The dashed line is the value for 
elementary pointlike bosons, A^ = 1/2. 



eq. (jl4p and we identify ii — a at the end of the computa- 
tion. An approximate solution of the functional integral 
for the effective action F gives access to thermodynamic 
quantities and correlation functions. 

The relation (fT4|) between a and n seems to be rather 
formal at this stage. In this section we will develop 
the physical interpretation of this formula. In this con- 
text the important distinction between microscopic and 
dressed molecules will appear. In a quantum mechanical 
computation for the physics of a Feshbach resonance the 
concept of dressed molecules arises from mixing effects 
between the open and closed channels. This is an impor- 
tant ingredient for the understanding of the crossover. 
Our functional integral formalism has to reproduce this 
channel mixing in vacuum and to extend it to the many 
body situation. In the functional integral formulation the 
quantum mechanical "mixing effects" are closely related 
to the wave function renormalization Z^. The concept of 
dressed molecules and their contribution to the density 
is directly related to the interpretation of eq. We 
stress however, that the functional integral evaluation of 
n for given cr = /i can be done completely independently 
of this interpretation. 



V. MOLECULE AND CONDENSATE 
FRACTION 

At this point the functional integral setting for the 
ultracold fermionic atom gases is fully specified. The pa- 
rameters flfl and can be extracted from a computation 
of two-atom scattering in the vacuum, taking the limit 
T ^ 0, n ^ [lB|. Rescaling with appropriate powers 
of kp yields the parameters c and h'^ for the many body 
system. The relation between a and n is determined by 



A. Exact expression for the bare molecule density 

The total density of atoms is composed of three com- 
ponents according to eq. ([7]), 

n = rip + '2nM + nc- (55) 

Here we have split up the contribution ng from ([7]) in 
a contribution from uncondensed bare molecules (con- 
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nected two-point function) 



flM = {(t)*4')~{(^*){<P)-hB 
= {(j)*(f))c- flB 



(56) 



and one from the condensate which only occurs below 
the critical temperature 



2hUP 



— ~ . 



(57) 



The ratio of atoms in the condensate arising from bare 
molecules is defined as 



the bosonic fluctuations can be interpreted as the den- 
sity of dressed molecules. Thus, again an estimate of the 
molecule density will be mandatory for the understand- 
ing of the crossover physics. A simple estimate of the 
dressed molecule density is the minimal ingredient be- 
yond mean field theory needed for a qualitatively correct 
description. We will refer to this as "extended mean field 
theory" . 

As a first step the evaluation of the molecule den- 
sity we may evaluate the classical approximation where 
the Yukawa interaction between open and closed channel 
atoms is neglected. This corresponds to the limit 

(for fixed a). Then G^?''\q^n) is the free propagator ^ 



n hi 



(58) 



~ "1/2 

and involves the Yukawa coupling ft,^ — 2Mh^/kp . 

The density of uncondensed bare molecules um ([55]) 
can be written in terms of the bare molecule propagator 
or connected two-point function 



nM{x)=T / dTG^{x,T\x,T) - fiB- 



(59) 



This formula is exact ^. It involves an additive renor- 
malization of a similar origin as for the fermionic 
density, cf. sect. IIII A| but with opposite sign due to 
Bose-statistics. In the homogeneous limit we can write 
eq. (|59p as a sum over integer Matsubara frequencies 



:G^{q,n) - ub- 



(60) 



The ratio of bare molecules to open channel atoms, 
fiM /np, is important for the understanding of the atom 
gas. Technically, it enters the field equation for the effec- 
tive chemical potential (j50p which involves npjn. In the 
regions in parameter space where um /n\s not very small 
a reliable computation of um is crucial for a quantitative 
understanding of the phase diagram. Such situations are 
e.g. realized in the BEC regime for narrow and inter- 
mediate resonances. However, for the crossover region 
for a broad Feshbach resonance (as for ^Li and ''"K) it 
turns out that um /n is small such that an approximation 
rip = n yields already a reasonable result. 

Nevertheless, the mean field approximation (jSip does 
not remain valid in all regions of the phase diagram. 
This is due to an additional cr-dependence of u aris- 
ing from the boson fluctuations which are neglected in 
MFT. We will argue below that this contribution from 



(61) 



Performing the Matsubara sum and inserting fiB = 
n/2 — 1/2 / (Pq/{2nY one obtains the familiar expres- 
sion involving the occupation numbers for bosons 



nM 



<Pq 
(27r)3 



AM 







exp 



2^JL. 



T 



(62) 



We note the role oi fiB for the removal of pieces that do 
not vanish for A — > cx). 



B. Fluctuation effects 

The fluctuation effects will change the form of G^{q, n). 
Details of the computation of G^ = V^^ are presented 
in app. [Cl In the symmetric phase we may account for 
the fluctuations by using 

P^iQ) = i2ninTZ^{a,T) + P^{q))5ab, 
P^iq) = A^ia,T,q^)q^ + rhlia,T). (63) 

We include here the fluctuations of the fermions (open 
channel atoms) for the computation of and Zrf,. The 
mass term will be evaluated by a gap equation which 
also includes the bosonic molecule fluctuations. In the 
superfluid phase {(f) ^ 0) the diagonalization of the prop- 
agator has to be performed carefully as discussed below 
in subsects. G,II. 

The gradient coefficient A^{a^T ^q^) is defined by eq. 
(Plj) and depends on the Yukawa coupling h^. For 
large it comes close to the classical value 1/4M. 
We plot the renormalized gradient coefficient A^icj^) = 
2MA^{q'^)/Zcf, for different c corresponding to the BCS, 



Eq. II59I I is valid for the normal phase T > Tc- For the super- 
fluid phase becomes a 2 X 2 matrix and the corresponding 
generalization of eq. II59I I is discussed in sect. IV HI 



* Note that we use 9 instead of Pa - this includes already the domi- 
nant fluctuation effects as motivated by the following paragraphs. 
Strictly speaking, the classical propagator features Pa- 
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FIG. 5: Weighted Bose distribution N^j^iq) = 
{exp{A^tf /T) — 1)~^ as a function of dimensionless momen- 
tum q = q/kp, for T — Tc- We show the three regimes; (a) 
BEC (c-^ = 1.5,1; = 0.218), (b) Crossover (c'^ = 0, fc = 
0.292) and (c) BCS (c"^ = -1.5, = 0.057). We com- 
pare our best estimate using A^{q) — A^{q)/Z^ (solid curve) 
with the approximation = A^{q = 0) (dashed curve) that 
we employ for the numerical estimates in this paper. Also 
indicated is the result for the classical gradient coefficient 
A^ = 1/2 (dashed-dotted) that strongly overestimates the 
molecule number in the crossover and BCS regimes. 



BEC and crossover regimes in fig. [H It is obvious that 
this gradient coefficient plays a major role. Large 
leads to an additional suppression of the occupation num- 
ber for modes with high g^, as anticipated in sect. IIII CI 

In the symmetric phase {(f)o = 0) the "mass term" is 



92 1/ 



rrij 



d(l)*d(j) 



4=0 



(64) 



We use here a language familiar in quantum field theory 
since the molecule wave (j) behaves like a massive field for 
TO^ > 0. (The propagator has a "gap".) For 77i| > 
the symmetric solution of the field equation, = 0, is 
stable whereas it becomes unstable for negative rfi^. For 
a second order phase transition the critical temperature 
therefore corresponds precisely to a vanishing mass term, 
m^{T = Tc) = 0. The mass term reads in the mean field 
approximation 



m 



{F)2 



d^UMFT 



2n 



(F) 



0=0 



(65) 



■AF)2 



For T = 0, (T = the definition of P d^H) implies : 
D — 2^. We note that the fermion fluctuations lower 
as compared to the microscopic term = — 2/i, cf. 
eq. (|19p . This effect enhances the occupation number 
for molecules - using ffi^ instead of would yield 
a much too small density of molecules! At the critical 
temperature the mass term vanishes (see below) . For 
T — Tc the fluctuation effects therefore concern the size 
and shape of ^0(9). 

As we have seen in our mean field computation in sect. 
mil the quantities Z^^A^j, and depend strongly on a. 
We may imagine to integrate first the fermion fluctua- 
tions in the functional integral ([6]) (with 77^ = rj — Q). 
The result is an intermediate "mean fleld action" for 
the remaining functional integral over bosonic fluctua- 
tions. The quadratic part (f>* {Q)(f){Q) in this action 
will be of the type of eq. ((63)) . Performing now the bo- 
son fluctuations will induce an additional contribution to 
dU/da = —rip. This will be interpreted below as the 
density of dressed molecules. 



C. Dressed Molecules 

Dressed molecules are quasi-particles with atom num- 
ber two. They are described by renormalized scalar fields 
<f> with a standard non-relativistic r-derivative in the ef- 
fective action. This allows for a standard association of 
the number density of quasi-particles with the correla- 
tion function for renormalized fields. With the effective 
action (fT3|) the relation between the fields for dressed and 
bare molecules reads 



.1/2; 



(66) 



Correspondingly, the dressed molecule density um be- 
comes 



riM = Z^riM, 



2nM 



(67) 
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and we find for large Zrf, a very substantial enhancement 
as compared to the bare molecule density um- 
We may define a renormalized gradient coefficient 



and mass term 



by 



A - 
^<pM — -77- 



or, in dimcnsionlcss units, 



A, 



Zrl, 



Za, 



2M 



TO. 



(68) 



(69) 



Then the quadratic part in the effective action for the 
bosons can be written in terms of (j) and m^^Atj,, with- 
out explicit reference to Zrp. (Similar rescalings can be 
made for any quantity entering our calculations. For a 
complete list of dimensionful, dimcnsionlcss and dimen- 
sionless renormalized quantities, cf. app. [F]) Using eq. 
in eq. (|60p the wave function Z^ can be factored out 



in Pij, such that 



Z^um 



(27r)3 



[exp ( 



-P0,fl(g) ' 
T ' 



A' 



(70) 



This is a standard bosonic particle number without the 
appearance of Z^, now expressed in terms of the effective 
renormalized inverse boson propagator (dimensionful and 
dimcnsionlcss version displayed) 



2MP^{q) 
kl Z^ 



= A4,W 



- A^f 



(71) 



The "dressed molecules" [3l|, l38j include both bare and 
fluctuation induced, effective molecules (cf. eq. (|38p ). 
We will see how the dressed molecule density um emerges 
naturally in the equation of state for the particle density 
below. 

We plot in fig. [5] the mode occupation number for the 
dressed molecules {q — \q\/kp) 



Nuiq) = (exp(P^,«(|g|)/T)-l)-i 
= (exp(P^(g)/f ) - l)-i 



(72) 



(weighted by a volume factor Sp). There we compare the 
"classical case" = /AM with the result including the 
fluctuation corrections. The normalization in the figure 
reflects directly the relative contribution to n^v/ 



1.3 

hp 

2^ 



d{\nq) fNM{q). 



(73) 



We observe a large fluctuation effect in the BCS regime 
(beyond the renormalization of to^). In this regime, how- 
ever, the overall role of the molecules is subdominant. On 




0.8 
0.6 
0.4 
0.2 




T = Tc 



(b) 



FIG. 6: Crossover at the critical temperature: Contribu- 
tions to the total particle number, showing the crossover 
from fermion to boson dominated physics, (a) Fractions of 
"dressed" densities in the large limit. We compare the 
results for two versions of the gap equation as in fig. [1] (b) 
Fractions of "bare" densities in the exact narrow resonance 
limit — > 0. Though the pictures are similar, the physical 
interpretation of the two plots differs as described in the text. 



the other hand, in the BEG regime the molecule distri- 
bution is rather insensitive to the details of the treat- 
ment of the fluctuations. The most important uncer- 
tainty from the "molecule" fluctuations therefore con- 
cerns the crossover regime. In the EEC and crossover 
regime the replacement of A^{q'^) by A^{q'^ = 0) results 
only in a moderate error. For simplicity we neglect the 
momentum dependence of ^^(q^) for the numerical re- 
sults in this work. 



D. Contribution of molecule fluctuations to the 
effective potential 

The computation of n^/ evaluates a one loop inte- 
gral which involves the molecule fluctuations. (Graph- 
ically, eqs. ([SO)) . ([5^ correspond to a closed loop for the 
molecule fluctuations with an insertion of a fi - deriva- 
tive.) A self-consistent approximation should therefore 
also include the effects of the molecule fluctuations in 
the computation of U and therefore (j)Q or Tc- Our func- 
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tional integral approach makes this necessity particularly 
apparent: The computation of the partition function 
Z involves the fluctuations of open channel atoms and 
molecules in a completely symmetric way (the variables 
'0 and (/) in eq. H])). There is no reason why the fluc- 
tuations of the fermionic atoms should be included and 
not the ones for the bosonic molecules. In particular, the 
critical region very close to Tc will be dominated by the 
boson fluctuations. 

For the effective potential the incorporation of the 
molecule fluctuations is achieved by adding to U the one 
loop contribution U\ from the fluctuations of </> 



U = Umft + Vf' =(y- '^IJ^W^ + V\"'+ Ul"'. (74) 

We can construct the bosonic contribution as the lead- 
ing order correction to the mean field result which omits 
boson fluctuations completely. For this purpose, we note 
that the fermion field appears only at quadratic order in 
the classical action ([T]). We can therefore integrate them 
out, which turns cq. ([5]) into a purely bosonic functional 
integral. 



(F) . jj{B) 



T[il; = 0, = - In j VS<f)exp ( - Smft[4' + 5(t>] (75) 

+j*S(f) + h.c.) 

with an intermediate action Smft depending on the field 
(j) ~ (j) + S(j). This is given by the exact expression 



1 



SmftI^] = ^;''l(/.]--lndet5W"^)[0] (76) 
= 4''^[0]-^Trln5«"^)[<^] 

where denotes the second variation w.r.t. the 

fermion fields. In the "classical approximation" one has 
r[0] = Smft[4>], while the one loop approximation corre- 
sponds to an expansion of Smft[4' + Scj)] to second order 



The mean field effective potential (|19p is obtained from 
eq. (|75p in the classical approximation. The next order 
contribution takes the Gaussian approximation for the 
fluctuations of the molecule field into account. In princi- 
ple, this requires the evaluation of highly nonlocal objects 
- the one-loop fermion fluctuations encoded in the Tr In- 
term in eq. ([75]) feature a complex frequency and mo- 
mentum dependence. However, since we are interested 
in the observable low energy properties of the system, 
we may apply a derivative expansion which only keeps 
the leading order terms in the frequency and momentum 
dependence. This precisely generates the wave function 
renormalization (|38l39p and the gradient coefRcient 
(|41IC13p . In this approximation, we find (up to an irrel- 
evant infinite constant and evaluated in the symmetric 
phase) the one loop result 



T 



(27r)3 



In 1 



'P<t>.B./T\ 



(77) 



where the spacelike part of the boson propagator is pre- 
cisely given by eq. (|7ip in the symmetric phase (for the 
result in the symmetry broken phase, cf. eq. (jlOSp ). This 
one loop formula has shortcomings and we will improve 
on it by solving appropriate gap equations in sect. IVIIII 
Nevertheless, it already contains the essential informa- 
tion for the different contributions to the density. 

Let us compare the effect of the /i- and cr-derivatives 
on the bosonic part of the effective potential, eq. (|77p . 
With the "classical" inverse boson propagator = 

/AM + V ~ 2^ one has dP^/dfj, = — 2 and we note the 
simple relation 



= -2nM. 



(78) 



On the other hand, the fermion loop corrections induce a 
a - dependence (at fixed fj.) of U{ , which contributes to 
fip. This contribution can be interpreted as the number 
of open channel atoms that are bound in the dressed 
molecules 



2nFM 



(B) 



da 



(79) 



The total number of dressed molecules is then given by 



(80) 



This could be taken as a possible alternative definition 

of 



(81) 



In the limit where the a - dependence of Zrp and Arj, can 
be neglected and dP^/d^ = — 2 one finds from 



(cf. eq. 



Za-1 = 



P^ A^q^ + 
) the relation 
dP^/da _ 1 dfhl 



(82) 



1 d^U 



dP^/dti 



2 da 



2 dadc/}* 



(83) 



We recover the MFT result in eq. ((40|) . The combined 
effect of the derivatives with respect to a and /x yields 
directly the number of dressed molecules 



riM 



7- ^fdu[ 



(B) 
1 



(B) 



da 



(84) 



The quantitative results shown in the figures are obtained 
by including the a - and (j) - dependence of P^ . They will 
be discussed in more detail in the next sections. 
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E. Open and closed channel atoms 

Some characteristic properties of the uhracold gas in- 
volve the fractions of open channel atoms, uncondensed 
bare molecules and condensed bare molecules 



ftp — ——, ^Ar — 



2n 



M 



n 



"c — — , 
n 



M 



1. 



(85) 



For example, the sum CIb = + measures the 
fraction of closed channel atoms, as observed in Q. The 
formal use of two different effective chemical potentials 
in the action 5*^ ([1]), i.e. a for the fermions and /i for the 
bosons, allows the simple association 



dU 



(86) 



and we recall that Qp also receives a contribution from 
u[^\ We have computed in [l^ and the results agree 
well with [s*] over several orders of magnitude. 

In the symmetric phase for T > Tc one has Qp + Qm = 
1. For small < 1 the BEC-BCS crossover is the 
crossover from small to large Ctp. In fig. [H] (b) we plot 
Clp and D,M as a function of the inverse concentration 
for T = Tc- The modifications of Z^, and fh?^ 
in the inverse propagator (|63p depend on the Yukawa 
coupling However, this influences the precise shape of 
the crossover between the BEC and BCS regime only for 
moderate >^ \. For smaller h^, the density fractions 
are insensitive to the fluctuation modifications. 

For the large values of encountered for the broad 
Feshbach resonances in ^Li and ""^K the contributions 
from the closed channel molecules VIm,^c become very 
small (cf. fig. [7](b)). The dressed molecules differ sub- 
stantially from the bare molecules (large Z^) and the 
crossover physics is better described in terms of dressed 
molecules. We display in fig. [5] (a) the fraction of dressed 
unbound atoms ^p and dressed molecules iljv/ for large 
values of and T = Tc- The fractions are not sensitive 
to the precise value of /i0 in the broad resonance limit 
hd, ^ oo, similar to the behavior found for small hd,. 



In fig. [7] (a) we plot the condensate fraction as 



function of 



for T = 0. We also show the fraction of 



closed channel atoms ^ in the condensate, flc, in fig. [7] 
(b). Both correspond to the broad resonance limit and 
we choose the Yukawa coupling appropriate for ^Li, ft,^ = 
610. 

For large one has fie ^ - indeed, the prob- 
ability Z^^ that a condensed di-atom state contains a 
bare molecule is small. In this case the major part of 
the condensate is due to open channel atoms, i.e. nc is 
dominated by a contribution from rip 



npc 



nc - nc ^ {Z^ ~ l)nc 



(89) 



nc + fi 



1 



-)np. 



Here the second line defines implicitly ftp hy the require- 
ment 



np + riM + nc = i. 



(90) 



(From nc + nM < n, flc + ^ <^ 1 one con- 
cludes Vtp 1^ and for low enough T one further expects 
CLm < ^c-) In the BCS limit this result is not surprising 
since we could have chosen a formulation without explicit 
molecule fields such that all atoms are described hy rip. 
(If the chemical potential multiplies only the total 
condensate fraction must be ilc = 1 ~ Hp ^°.) 

The total number of open channel atoms can therefore 
be found in three channels, rip = np+npc + 2npM- Here 
np denotes the unbound dressed fermionic atoms, 2npM 
the open channel atoms contained in dressed molecules 
and npc the ones in the condensate. As long as the 
fermionic and bosonic contributions to U can be sepa- 
rated we have the identities 



np,o 



nM 



np 



npc = 



1/dUl 



da ' 



(91) 



da 



and 



F. Condensate fraction 

The total number of atoms in the condensate depends 
on the expectation value of the renormalized field <f>ji = 

1 /2 ~ 

— Zj 0. We will mainly use the dimensionless field 



,-3/2,7 , ,-3/2„l/2- 
= kp Wr) = kp ZJ C, 



(87) 



With p = 4>* (j), pq = 4>q4>o the dressed condensate fraction 
can be defined as 



Vie = = 2kp4>Q4)Q/n = 67r po = Z^Vlc- 



(88) 



riF.o + 271m + nc 
np + 2nM + nc 
lip + 2nM + nc- 



(92) 



The definition of a condensate fraction in terms of the 
superfluid order parameter po is rather simple and ap- 
pealing. Nevertheless, this may not correspond precisely 



® The fig. for includes renormalization effects for discussed 
in [11. 

The difference between Qc a-nd 1 — Q,p can be traced back to 
the appearance of the chemical potential ^ in the effective four- 
fermion interaction mil l. 
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to the condensate fraction as defined by a given exper- 
imental setup. The ambiguity is even larger when we 
come to the concepts of uncondensed fermionic atoms 
and molecules. The distinction becomes somewhat arbi- 
trary if we include higher loops for the computation of it, 
where bosonic and fermionic fluctuations arc mixed. As 
an example for the ambiguities in the definition one may 
try to extract the number of open channel atoms in the 
condensate directly from the (j) - dependence of np. We 
can decompose the fermion contribution into a part for 
vanishing condensate and a "condensate contribution" 
due to (f> ^ 



(F) 



The association 



(F) 



npc 



dAUl 



(F) 



da 



yields an alternative definition of Z^, 

1 9A[/}^' 



Z'-l 



flc 



24)* 4) da 



(93) 



(94) 



(95) 



With this definition the number of unbound atoms np 
becomes 



np — np ~ npc — npM 



du[''\4> = o) 



which amounts to the standard number density in a 
Fermi gas. We note that Z'^ ((95| coincides with (ISS)) 
in the limit where d^U/dad(f>*d(f> is dominated by the 
term quadratic in (j) in U[^\ From fig. [7] (a) we see that 
this is the case in the BEC and BCS regimes. 



G. Excitations in the superfluid phase 

The bosonic excitations in the superfluid phase are 
analogous to a purely bosonic theory as for superfluid- 
ity in ^He [sij . Due to the nonvanishing order parameter 
the matrix for the inverse renormalized propagator 
contains now "off-diagonal" entries form terms ~ Po4'4' 
or pq4>*4>* where po — kp^Z^{(j))*{(j)). It is convenient to 
use a basis of real fields (t>i , (t>2 



(97) 



0= -^(01 +102), P=\^ 

In this basis and = V^^ are 2x2 matrices and the 
exact eq. ((60)) for the bare molecule density is replaced 
by 

T 



nM{x) = — tr / dTG^{x,T;x,T) 



riB- 



(98) 



We expand in the superfluid phase around the minimum 
of the potential at p ~ pq 



A. 



ip-Po? 



(99) 
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FIG. 7: (a) Contributions to the total particle density for 
T = in the large limit: The fraction of dressed molecules 
flu (dashed line) is largest in the crossover regime. The con- 
densate fraction Qc (solid line) grows to one in the BEC 
regime. The solid line corresponds to flc ~ Z^Q.c whereas 
the dashed-dotted line uses Z'^ (eq. ((95)l) instead of Z^. (b) 
Fractions of the bare or closed channel molecules. In contrast 
to the dressed molecules, they are 0{h~^^). The dominant con- 
tribution arises from the condensed bare molecules Q,c (solid 
line). The contribution from noncondensed bare molecules 
Q.M at T = remains very small. 



such that the mass matrix 



' d4ad<f>b 



P=Po 



becomes 



V 0, 



(100) 



(101) 



Without loss of generality we have taken here the order 
parameter in the 0i direction, ^i.o — V^Po, 4'2.o — 0, and 
we recognize the flat direction in the potential (vanishing 
eigenvalue of in the "Goldstone direction" 4>2- In 
contrast, the "radial mode" (f)i has a nonvanishing mass 
term 2X^pQ. 

In the basis (0i, 42) the term containing the r - deriva- 
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tive is off-diagonal (neglecting total derivatives) 

J dT(j)*dr4> = i j dTcPxdrh ^ dT(t)2dr4>l. (102) 

In momentum space one therefore finds for the renormal- 
ized inverse propagator = 2M / {Z^kp)V^'^^ 



Co , A^q^ 



(103) 



where we use eq. (j69l) and the renormalized order pa- 
rameter and quartic coupling 



(104) 



(For a list of the relations between dimensionless and di- 
mensionless renormalized parameters cf. app. [F]) This 
has an important consequence: The propagating exci- 
tations correspond to frequencies w which obtain from 
the Matsubara frequencies t^s by analytic continuation 
u)B iw. This corresponds to the analytic continua- 
tion from Euclidean time r to real or "Minkowski" time 
t = —IT. Now has a pole or a zero eigenvalue. 
The eigenvalues A of therefore obey 

{A^f + 2X^po ~ X){A^f ~X)-Cj^^O. (105) 

Vanishing eigenvalues A therefore lead to the dispersion 
relation 



(D^ = A^f{2X^po + A^f). 



(106) 



For small this yields the linear dispersion relation char- 
acteristic for superfluidity 



Cu = ^2A^X^poVf: (107) 
from which we can read off the speed of sound 

(108) 



Vs = fcF/(2M) y/2A^X^po. 



H. Molecule density in the superfluid phase 

The density of dressed molecules in the superfluid 
phase obeys 



T r (pq 



n 



d^q A^q"^ + X^pq 

{2TT)^u,l + A^q^{A^q^ + 2X4 



(109) 



where the Matsubara summation can be performed ana- 
lytically 



A^q^ + X^po 



fcf f d^q f 

2 J (2^)n^A^q2(^^q2 + 2A^po) 

xcoth^5!W±!5M_i|. (110) 
2T i ^ ' 



Due to the subtraction of (the term — 1 in the curled 
bracket) the momentum integral is UV finite. 

We have used dimensionless units and may introduce 

_ J {A^q'^ + m'^) / {2T) symmetric phase 



\ X^po/(2f) 
where (for SSB) 



superfluid phase 

symmetric phase 
superfluid phase 



(111) 



a0 = ^A^q^{A^q^+2X^po)/{2f) (112) 

— \/ + 2Ka 

is the bosonic analog to eq. ([^0)1 . This allows us to write 
the dressed molecule density in both phases as 



(SYM) _ 
{SSB) _ 



(2 



n 



M 



— -:7 expia 
d^q 



1 



{2TTf 



/ a + K , \ , 
cothofrf,- 1 . 113) 



Note the structure r2 ^ 1/2 f 6'^ {-Q)V4,{Q)<f){Q). 



At the phase boundary the two definitions coincide since 
= a {rn^ = 0, po = and k — 0). 

The equations (|113p for the density of dressed 
molecules involve the renormalized coupling A^. This 
describes the full vertex and therefore A^ is a momen- 
tum dependent function. This also holds for Ac/,. We 
will neglect the momentum dependence of A^, as moti- 
vated by fig. [5] for T — Tc- For A0 this issue is more 
involved since X^{q) vanishes for q —f Q {T < Tc) due 
to the molecule fluctuations, as shown in app. [E] We 
observe that for small T the momentum integration in 
eq. piOp is dominated by the range « X^pa/A^. The 
infrared suppression of X^{q — > 0) is therefore not effec- 
tive. For the density equation we will simply omit the 
contribution of the molecule fluctuations to A^ and ap- 
proximate A0 ~ X^^^ with A^^^ evaluated at q^ — 0. 

The inclusion of the molecule density um is impor- 
tant for quantitative accuracy even at T = 0. This is 
demonstrated in fig. [5] where we show the crossover for 
the effective chemical potential a and the gap A as a 
function of . The agreement with quantum Monte 
Carlo simulations [IJ is substantially improved as com- 
pared to mean field theory, and also compares reason ably 
well with other analytical approaches (Strinati et al [52l |. 
A = 0.53, a = 0.445). 
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FIG. 8: Solutions of the coupled gap and density equations 
at T = 0, (a) efltective chemical potential, (b) gap parameter 
A — \ff. We also show (dashed)the mean field result (ob- 
tained by setting the bosonic contribution nu = in the den- 
sity equation). Our result can be compared to QMC calcula- 
tions [37| performed at c"^ = which find a = 0.44, A = 0.54. 
Our solution yields <j = 0.50, A = 0.53, and improves as com- 
pared to the MFT resuh a = 0.63, A = 0.65. 



VI. EFFECTIVE FIELD FOR ATOM DENSITY 



Before proceeding in sect. I Villi to a description of our 
computation of the effects from the molecule fluctuations 
on the effective potential we present in this section an 
improvement of our formalism which treats the fermionic 
and bosonic fluctuations in an even more symmetric way. 
Indeed, in the thermodynamic equihbrium the molecule 
propagator should involve the same effective chemical po- 
tential a as the propagator of the unbound atoms. (So 
far it involves /i instead of tr.) 

If one is interested in the separate contributions from 
open and closed channel atoms one may formally con- 
sider different chemical potentials tr and /i multiplying 
and in the action. Then np and fiB can be 
associated to the variation of the effective action with re- 
spect to a and /i. At the end of the computations one 
has to identify a — [i). For many purposes, however, one 
only needs the total number of atoms n. It seems then 
advantageous to modify our formulation such that the 



field cr is associated to n instead of np. In this section, 
we treat cr as a classical field such that its role is reduced 
to a source term a = J . In app. |^ we will consider 
the more general case where a is treated as a fluctuating 
field. 

Identifying <j — ji in the bare action ([T]), the cr- 
derivative of the effective potential generates the total 
particle number density, 



9^ 



np + 2n 



M 



20* 



(114) 



In the presence of interactions it is hard to evaluate the 
full bare correlation functions in the above form explic- 
itly - it premises the solution of the full quantum field 
theory. As we have discussed in sects. IV D[ IV El it is 
more practicable to decompose n according to 



MFT ■ 



c/(-)) 



da 



nc + npfi + 2nM (115) 



which is a "mixed" representation involving both bare 
and renormalized quantities. In the next section we will 
see that this equation reduces to an equation of state for 
"fundamental" bosons in the BEC regime, which involves 
dressed quantities only (eq. (|12ip ). 

Beyond the classical bosonic propagator we have im- 
plemented the cr - dependence of Ui^'^ by the a - de- 
pendence of mi. This yields 



1 

2~d^ 



1 d^U 



2 d(jd(j)*d(l) 



(116) 



where we note that the contribution from the bare 
molecules is already included in the classical a - depen- 



dence of , leading to the replacement 
eq. ()83p . In practice, we use the approximation 



1 in 



^_i d^uf^ 



2 dac 



(117) 



This has a simple interpretation: due to the fermionic 
fluctuations the classical "cubic coupling" —2a(t>*(j) is re- 
placed by a renormalized coupling —2Z^a4>*4>. We can 
compute UM directly from eqs. (|73l) . (ffS]) . (TfOl) . ([63 ]) . 



VII. BEC LIMIT 



In sect. IV Gl we have deflned a multiplicative renormal- 
ization scheme for the bosons by rescaling all couplings 



A more rigorous derivation of the equation of state via the 
Noether construction for the conserved charge associated to the 
global (7(l)-symmetry reveals that this choice is uniquely fixed 
by requiring consistency within our truncation (linear frequency 
dependence) O . 
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FIG. 9: Gradient coefficient = Atp/Z^ for T = 0, as a 
function of the inverse concentration. In the BEG limit, A^f, 
takes the classical value 1/2 as appropriate for elementary 
bosons of mass 2M. 



and fields with an appropriate power of such that the 
timelike derivative in the effective action has a unit coef- 
ficient. This leads to the notion of dressed molecules in 
our formalism. In this section we discuss the implications 
of this prescription in the BEC regime. Here we will see 
that the dressed molecules behave just as "elementary" 
bosons. 

The propagator for the dressed molecules has still a 
nontrivial renormalization factor for its dependence on 
the spacelike momenta (|63)) . The ratio = A^jZ^ is 
shown as a function of c in fig. [51 It is instructive to eval- 
uate A^ in the BEC limit where the propagator should 
describe the propagation of dressed molecules. The in- 
tegrals in A^ and Z^ can be evaluated analytically (cf. 
app. ini eqs. (|D3| ID4|) ) irrespective to the thermody- 
namic phase of the system. The universal result does 
not depend on /i^, 



A. 







(118) 



This makes it particularly clear that for hoth the limit 
of small /i0 — > and the limit of large we re- 
cover composite, but pointlike bosons with wave func- 
tion A^ = 1/2, corresponding to a kinetic energy f? j^M. 
However, in the first case of small we deal with micro- 
scopic closed channel molecules, while in the second case 
of large the emerging bosons are open channel pairs, 
however behaving just as if they were pointlike bosons in 
the BEC limit. For large we could equally well drop 
the classical piece in A^ and Z^. This corresponds to a 
purely pointlike interaction for the fermions without any 
explicit reference to molecules at all, cf. eq. (fTTT) . and 
the discussion in [l^. Nevertheless, bound atom pairs 
emerge that behave just as pointlike particles. 

A similar aspect of universality is found for the four- 
boson coupling = A^/Z|. In the approximation where 
we neglect bosonic contributions to the four-boson cou- 
pling and to the gap equation for the mass (BCS gap 



equation), we find = = Stt/v— ct such that the 
bosonic scattering length becomes au = 2a/j. This is 
the Born approximation for the scattering of the com- 
posite bosons. However, more accurate treatments on 
this issue have shown aujo-R — 0.6. This result is ob- 
tained from the solution of the four-body Schodinger 
equation fssj and form a numerically demanding dia- 
grammatic approach l5j|, and has also been confirmed 
in QMC simulations [53. A resummation of the effec- 
tive boson interaction vertices in vacuum as done in fHB] 
yields aula-R = 0.75(4). Our present Schwinger-Dyson 
approach does not correctly account for the real situation 
since the contribution from molecule fiuctuations in the 
vacuum is missing. In the frame of functional renormal- 
ization group equations, this deficiency has been reme- 
died. Simple truncations yield qm/cir = 0.71 — 0.92 [49|. 
It is instructive to investigate eq. (|113p in the limit 

c 

2aR. 

then coincides precisely with the Bogoliubov formula for 
bosons of mass 2M and the above scattering length. 



^ — > 00 (BEC limit) where we have 1/2, ~ 

The density in the superfiuid phase um (|113p 



(2^ 



exp 2q!0 — 1 



(119) 



where we used the relations connecting the Bogoliubov 
transformation coefficients and our expressions (jllll 

nil. 



1 ,a + K 
2' 



1), + 



(120) 



We emphasize again that we can perform first the limit 
00 where we recover a purely fermionic model with 
pointlike interaction and no explicit molecule degrees of 
freedom (cf. [3). Subsequently wc may consider large 
where the approximations leading to eq. (|119p be- 
come valid. This shows that the Bogoliubov formula for 
weakly interacting "fundamental" bosons can be recov- 
ered from a purely fermionic model! In our approach, this 
result emerges in the simultaneous limit c^^ — > 00 (BEC 
regime) , /i^ — )■ 00 (broad resonance regime). A similar re- 
sult has been established by Strinati et al. p!8l.[l9l.[20l.l5^. 
who work in a purely fermionic setting, or, in our lan- 
guage, in the broad resonance limit ^ 00 from the 
outset. 

This observation is strengthened by an investigation of 
the dimensionless density equation (|115p at T = in the 
BEC hmit: 



1 = 37r2(- 



37r2(2Z^^ 



2nM) 



2nM) 



M- 



(121) 



Here we have dropped the term fie = 0{h^'^) and we 
may similarly neglect ^m- In the first line of eq. (|12ip 
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the first term is the explicit result for Qp.o- The sec- 
ond line uses the explicit result for (|D4p and we re- 
cover the definition of flc (|88p . This is precisely the 
density equation one obtains when assuming fundamen- 
tal bosons of mass 2AI. This is an important result: 
Our Z,f, - renormalization procedure generates precisely 
the macrophysics we would have obtained when start- 
ing microscopically with a purely bosonic action. In our 
approach however, the bosons emerge dynamically. 
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VIII. GAP EQUATION FOR THE MOLECULE 
PROPAGATOR 

The formulation of the problem as an Euclidean func- 
tional integral for fcrmionic and bosonic fields permits 
the use of many of the highly developed methods of quan- 
tum field theory and statistical physics. It is an ideal 
starting point for systematic improvements beyond MFT. 
We have mentioned in the previous sections that one pos- 
sible alternative to the standard one loop approximation 
could be to first integrate out the fermions and then per- 
form the remaining (/i-integral in one loop order. This 
procedure has, however, some problems. 

Let us consider the dimensionless renormalized inverse 
bosonic propagator (j63|) {ub — 0) in the approximation 



The exact mass matrix 



'•0,ofc 



0=0 



(122) 



(123) 



is diagonal but in general non-degenerate in the 0i , 4>2 
basis. It vanishes at the critical temperature of a sec- 
ond order phase transition. Already after the first step 
of the ■0 - integration Vcf, differs from the classical in- 
verse molecule propagator. The molecule propagator will 
appear in the molecule fluctuations as described by the 
partition function 



Z 



V(t>e- 



(124) 



where Smft is the intermediate action resulting from the 
Gaussian integration of the fermion fields. It is composed 
of the classical boson piece and a loop contribution from 
the fermion fields, 



S 



MFT 



= d^x{4)*{Z^[4>\dr~A^[4)\^ + v-2a)^ 



iF)l 



+ ...}. 



(125) 



At this stage the contribution from the fermion loop still 
depends on the fluctuating boson field. For example, the 
formula for the dressed molecule density in the symmetric 
phase 



(2^) 



7r)3 V 



P<pig)/T 



1 



IdUl 



(B) 



2 da 



(126) 



FIG. 10: Lowest order Schwinger-Dyson equation for the in- 
verse molecule propagator (double dashed line) in the sym- 
metric phase. The first two terms on the rhs denote the "mean 
field inverse propagator" after integrating out the fermionic 
fluctuations with a dashed line for the classical inverse 
molecule propagator and a solid line for the fermion prop- 
agator. The third term on the rhs accounts for the molecule 
fluctuations. Here is the molecule self-interaction in- 

duced by the fermion fluctuations. 



involves (|122[) . It generalizes eq. (|62p with fj, re- 
placed by a. We see that the inverse propagator ap- 

pears in the bosonic fluctuation contribution [/{ to the 
effective potential and influences the field equation for 0. 

For a simple one loop evaluation of the functional in- 
tegral (|124p one would have to use a propagator with 

TO^^-* ^ instead of m| in eq. (|122p . This has an important 
shortcoming. Due to the difference between Umft and U 
the masslike term vanishes at some temperature 

different from T^. As a consequence of this mismatch 
one observes a first order phase transition for sufficiently 
large values of h^. Clearly such a first order phase tran- 
sition may be suspected to arise from an insufficient ap- 
proximation. (A similar fake first order transition has 
been observed for relativistic scalar theories. There it 
is well understood [stI [58| that an appropriate resum- 
mation (e.g. by renormalization group methods) cures 
the disease and the true phase transition is second or- 
der. For the crossover problem the second order nature 
of the phase transition has recently been established by 
functional renormalization group methods for the whole 
range of concentrations, cf. [ioj ) 

In order to improve this situation we use for the 
bosonic fluctuations Schwinger-Dyson type equations 
where the inverse propagator P^ involves the second 
derivative of the full effective potential U (rather than 
Umft)- It is obvious that this is needed for a reliable esti- 
mate of um since the exact expression (1591) and therefore 
also (|126p involves the full molecule propagator includ- 



Eq. H126P approximates the exact eq. H60| l in the limit where cor- 
rections to the dependence of on the Matsubara frequencies 
can be neglected, see app. [C] The second part can be viewed as 
a Schwinger-Dyson equation for the a - dependence of U. 



23 



ing the contribution of the molecule fluctuations. Since 
our functional integral treats the molecule fluctuations 
exactly on the same footing as the fermionic atom fluc- 
tuations we can derive the lowest order Schwinger-Dyson 
or gap equation for in the standard way. For the 
symmetric phase it is graphically represented in fig. [10] 
and involves the full propagator and therefore 



A. Symmetric phase 

We will consider a Taylor expansion of u = 2Mkp^U 
in terms of the invariant 



(127) 



Correspondingly we use a renormalized Yukawa coupling 
and four-boson vertex, 

= 2Mh^/{kFZ^f/^, = X^/Z'l. (128) 

For the symmetric phase we expand 



(129) 



1 



UMFT = 'm^p+ -\\ ' p 



The gap equation for takes the form 



2 {F)2 ^ (B)2 



(130) 



AB)2 



Here A^ is the effective vertex involving four bosonic 
fields ^ (0*0)^, as induced by the fermion fiuctuations. 
It depends on and Neglecting the Wn - dependence, 

i.e. replacing xf \q,Cjn) ^ >^f\q) = >^f\q,'^n = 0), 
one can perform the Matsubara sum 



(131) 



AB)2 



(27r)3 > 



2f 



(132) 



We have not yet computed the momentum dependence 
of A^ (q) but a simple qualitative consideration of the 
relevant diagram shows that for large q^ one has a fast 
decay X''^\q) cx q~^. This makes the momentum inte- 
gral (|132p ultraviolet finite. It will be dominated by small 
values of q^. For our purpose we consider a crude approx- 

( F) ( F) { F) 

imation where we replace A^ (9) ^ •^^ (9 = 0). 

Of course, we have now to restrict the momentum inte- 
gration to low momenta. This can be done efficiently by 
subtracting the leading UV divergence similar as in the 



computation of Hm, i-e. replacing cothx by cothx — 1. 
This procedure yields 



AB)2 



2A 



(F) 



exp 2a 



1 



(133) 



We recognize on the rhs of eq. (|133p the expression for 
the number density of dressed molecules and obtain the 
gap equation 



(-F')2 , V- 



(134) 



where we recall that flM depends on to^. Our gap equa- 
tion has a simple interpretation: The bosonic contribu- 
tion to m| vanishes in the limit where only very few 
dressed molecules play a role {^Im 0) or for vanish- 
ing coupling. We are aware that our treatment of the 
suppression of the high momentum contributions is some- 
what crude. It accounts, however, for the relevant physics 
and a more reliable treatment would require a quite in- 

(F) 

volved computation of A^ {q,Lb). This complication is 
an inherent problem of gap equations which often require 
the knowledge of effective couplings over a large momen- 
tum range. As an alternative method one may employ 
functional renormalization (soj for the crossover problem 
[4^, where only the knowledge of couplings in a narrow 
momentum interval is required at every renormalization 
step. 

( F) 2 

The fermionic contribution to the mass term rrr^ 
(cf. eq. ([S5)) 1 reads exphcitly 



m 



(F)2 



(F)2 



2ct du\ 



dp 
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Zrh 



/?2 

AT 



d^q 
(27r)3 



1 2T 

— tanh70 37 

i4> r ■ 



The expression in the last line has to be evaluated with 
70 = 7 in the symmetric phase. The r.h.s. of the 

gap equation (|134[) involves the coupling A^^' for the 
molecule- molecule interactions 



,F) _ 2MkF^(F) _ d^uf^ 
Aa — — — Ai — 



32f3 



9p2 



(136) 



d^q 



■^—^ {70 ^ tanh 14,-1/ cosh 7^ } . 



Again, in the symmetric phase A^ has to be evaluated 
at p = 0, i.e. 70 = 7. 

The molecule fiuctuations give a positive contribution 
to 7710, opposite to the fermionic fluctuations. This has 
a simple interpretation. The fermion fluctuations in- 

(F) 

duce a self- interaction between the molecules ~ A^ . In 
turn, the fluctuations of the molecules behave similarly 
to interacting fundamental bosons and modify the two 
point function for the molecules. The quantum correc- 
tions to the fermionic and bosonic fluctuations to the in- 
verse molecule propagator are represented graphically in 
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FIG. 



11; Temperature dependence of the gap A = \/F at the 
resonance. The role of molecule fluctuations and the uncer- 
tainties in their treatment for T —> Tc are demonstrated by 
four choices of in the gap and density equation. The crit- 
ical temperatures are indicated by vertical dashed lines, with 
values fc = 0.259, 0.292. 



fig. [TOl We emphasize that an additional microscopic 
molecule interaction could now easily be incorporated 



■'classi- 



by adding to the mean field value for A 

cal part" A^'^'"'. In this case the renormalization of Pa 
discussed in sect. IIII Al would be modified. For a con- 

(cl) 

stant A^ the UV - divergent part would contribute to 
miiT = 0,71 = 0). One would again end up with a con- 



tribution of the form (|133p . now with A^^-* replaced by 

^^^'"'^ ■ Our approximation (|133p therefore treats the 

interactions of dressed molecules similar to fundamental 
interacting bosons. 

Actually, the Schwingcr-Dyson equation for the 
molecule propagator (fig. [10]) also describes the contri- 
bution of molecule fluctuations to the momentum depen- 
dent part encoded in A^. In the limit of a momentum 
— t p\ 

independent A^ the contribution of the boson loop to 
vanishes in the symmetric phase. 



B. Superfluid phase 

In the superfluid phase we choose an expansion of u 
around the minimum at 



namely 



(137) 



(138) 



and define correspondingly 

duMFT 



AF)2 



dp 



d'^UMFT 



Pa 



Pa 



(139) 



Our truncation of u[^'^ approximates 



AB)2 



(p- po) + 2 A 



ip-Po? 



(140) 



The location of the minimum po is determined by the 
condition 



AF)2 



AB)2 



0, 



AB)2 



r.(B) 



dp 



P=Pa 



(141) 



This defines the gap equation for po, which is the equiv- 
alent of eq. (|134p for the superfluid phase. The com- 

( B\ 2 

putation of the bosonic contribution rri^ encounters 
the same problems as for rri'J^^ ^ in the symmetric phase. 



Again we replace X^^\q) by a constant A^'^'' evaluated 
for q = and subtract the leading UV divergence of the 
momentum integral. This results in 



,iF) 



AB)2 



2\\ 



(F) [ d^q A^q^ + A0PO/2 



(27r)3 ^A^q^{A^q^ + 2\^po) 



(142) 



x(expy^^0g2(A0g2 + 2A0Po)/T- 1) . 

(I111I111I112|) we ar- 



In terms of the shorthands a, k 
rive at the gap equation for po 



AF)2 



2X\ 



(F) f (Pq_a + K/2 



(27r)3 



exp 2a^ 



1)"' =0. 



(143) 



The quantity contains a mass term 2X^pQ which 



involves the "full" vertex Arf 



A 



{F) 



+ Xi^\ We have 



computed the Schwinger-Dyson eq. for A^ in app. [E] For 
zero momentum g ^ we find that A^ vanishes in the 
superfluid phase. This is, however, only part of the story 
since the gap equation involves a momentum dependent 
vertex X(j,{q). In order to demonstrate the uncertainty 
arising from our lack of knowledge of X^{q) we present in 
fig. [m our results for different choices of A^ in the gap 
eq. dug). 

In detail, we show in fig. [Tl]four approximation scenar- 
ios; (i) the "standard" BCS gap equation (long dashed) 

neglects the molecule fluctuations, i.e. A^ = A^^-* = 
in both the gap and the density equation. This yields a 
second order transition but disagrees with QMC results 
for r — > 0. (ii) Bogoliubov density (short dashed) with 

(F) 

A^ in the density equation, while the molecule fluctu- 
ations in the gap equation are neglected. This improves 
the behavior for T ^ 0, but induces a fake first order 
phase transition for T ^ Tc- (iii) Neglection of molecule 
fluctuations in the effective coupling A^ (dashed-dotted), 

(F) 

i.e. we use A^ = A^ in the density and gap equation, 
(iv) Our best estimate (solid line) includes also correc- 
tions from molecule fluctuations for A^. As described in 
this section we use A^ in the propagator of the diagram 
in the gap equation, whereas the coefficient multiplying 
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the diagram is given by A^^'' . (We use A^^^ in the density 
equation.) The first order nature is weaker than in (iii), 
but stiU present. In a recent functional renormalization 
group treatment we have established the second order 
nature of the phase transition throughout the crossover, 
indicating that we now control the universal long range 
physics governing the phase transition [49| . 



C. Phase transition 

For T ^ Tc the gap equations pM]) and pl5)) match 
since po = 0, = a, k — 0. Also the expression for the 
molecule density becomes particularly simple 

3r(3/2)c(3/2)/rx3/2 

= ^2 [j-J ■ (144) 



IX. CONCLUSIONS 

Our functional integral investigation of ultracold 
fermionic atoms strengthens the picture of a smooth 
crossover between Bose-Einstein condensation of 
molecules and BCS type superfluidity. One and the 
same field (j> can both describe molecules and collective 
atom excitations of the type of Cooper pairs. The two 
different pictures correspond to two regions in the space 
of microscopic couplings and external parameters (T, n) . 
In dependence on a concentration parameter c, which 
is related to the magnetic field B, one can continuously 
change from one to the other region. Away from the 
critical temperature the "macroscopic observables" 
are typically analytic in c^^ despite the divergence 
of the scattering length for two-atom scattering at 
= 0. Only for T ^ one expects to encounter the 
non-analytic critical behavior which is characteristic for 
a second order phase transition in the universality class 
of 0(2) - Heisenberg magnets. 

For small and moderate Yukawa couplings the phys- 
ical picture of the crossover for T > Tc\s rather simple. 
Far on the EEC side the low temperature physics is dom- 
inated by tightly bound molecules. As the concentration 
increases, the molecule waves start to mix with collective 
di-atom states. More precisely, the molecule state with 
momentum p mixes with pairs of atoms with momenta 
(fi , q2 that are correlated in momentum space such that 
9*1+92 = p (Cooper pairs). On the BCS side the relevant 
state becomes dominantly a Cooper pair. Our formalism 
uses only one field (j) for both the molecule and Cooper 
pair states. Nevertheless, the relative importance of the 
microscopic molecule versus the collective Cooper pair is 
reflected in the propagator of 0, i.e. the wave function 
renormalization Zip and the gradient coefficient A,p. 

On the BEC side the propagation of corresponds to 
free molecules and takes the classical value A^^''^ = 
1/4M, and similarly for the wave function renormaliza- 



tion, =1. The increasing mixing in the crossover re- 
gion results in and Z,p growing larger than A'^^^^ , z'^^^^ . 

Indeed, the diagram responsible for the increase ofAp, Z^f, 
precisely corresponds to a molecule changing to a virtual 
atom pair and back. Finally, far in the BCS regime one 
has A,p A''^'^ , Zcj, ^ Z^'^ such that the classical con- 
tribution is negligible. The existence of a microscopic 
molecule is not crucial anymore. At long distances the 
physics looses memory of the microscopic details - the ap- 
proximation of a pointlike effective interaction between 
atoms becomes valid. 

For a broad Feshbach resonance the situation is, in 
principle, similar. However, the region where Zp and 

A^/A^*^'^ are large covers now both sides of the Feshbach 
resonance. The difference between the bare (microscopic) 
and dressed molecules is a crucial feature. For a broad 
Feshbach resonance the microscopic molecules are irrel- 
evant for the whole crossover region. For this limit we 
have demonstrated in sect. IVIII that our formalism is 
equivalent to a purely fermionic model with a pointlike 
interaction. This demonstrates that in the broad reso- 
nance limit, the use of the two channel model is rather a 
question of computational ease than a physical issue - it 
describes the same physics as a single channel model. 

On the other hand, the two channel model is more 
general and also covers the case of narrow resonances. 
This issue is discussed in more detail in [T5| . It remains 
to be seen if it will be possible to investigate the narrow 
resonances also experimentally in the future. 

The universality of the low temperature physics and 
the crossover can be traced back to the "renormaliz- 
able couplings" of an effective non-relativistic quantum 
field theory for long distances. More precisely, the long 
distance physics will only depend on the relevant (or 
marginal) couplings in the infrared renormalization flow. 
These are precisely the dimensionless parameters c and 
hp. Improved approximations will influence the relation 
between microscopic atomic and molecular physics and 
{c,h,p). At this point also "subleading interactions" like 
the CT-exchange or other interactions contributing to the 
background scattering length abg will play a role. 

However, universality predicts that the relations be- 
tween long-distance ("macroscopic") quantities should 
become computable only in terms of the "renormaliz- 
able couplings" (or c) and /i^. In consequence, we 
find a universal phase diagram in terms of these parame- 
ters, where the "memory" of the microscopic physics only 
concerns the values of c and /i^, providing for a universal 
phase diagram in terms of these parameters. In addition, 
the universal critical exponents and amplitude ratios for 
T ^ Tc will be independent of c and hrp. Furthermore, 
the BCS limit is independent of /i^ since only one param- 
eter (c) characterizes the effective pointlike atom inter- 
action. Also the BEC limit does not depend on h^p since 
the fluctuations of unbound atoms become irrelevant. On 
the other hand, the behavior in the crossover region can 
depend on h,p as an important universal parameter. This 
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demonstrates that a pointlike approximation for the ef- 
fective atom interaction is not always applicable for small 
\c-'\. 

The couplings and /i^ are related by the multiplica- 

1/2 

five factor . In the broad resonance limit —> oo 
one finds that the renormalized coupling h(j, is given by 
an infrared fixed point. It therefore ceases to be an inde- 
pendent coupling. 

For a more general form of the microscopic action we 
expect that deviations from universality become most 
visible far in the BCS regime, where further channels 
may play a role for the effective interaction, as well as 
in the far BEC regime, where more details of the mi- 
croscopic dispersion relation for the molecules and their 
microscopic interactions may become relevant. Going 
beyond our particular ansatz for the microscopic action 
universality should actually work best in the crossover 
region of small \c~^\. On the other hand, for a given mi- 
croscopic action the results are most accurate in the BEC 
and BCS regimes where fluctuation effects are most easily 
controlled. It is precisely the presence of strong fluctu- 
ation or renormalization effects in the crossover regime 
that is responsible for the "loss of microscopic memory" 
and therefore for universality! 

The formulation as a Yukawa model solves the prob- 
lems with a large effective scattering length for the atoms 
in the crossover region in a simple way. The large scat- 
tering length is simply due to the small "mass term" 
of the molecule field. It does not require a large atom- 
molecule interaction oc h^. As long as the dimensionless 
Yukawa or Feshbach coupling remains small, the re- 
gion of small and vanishing \c~^\ (large or diverging ef- 
fective atom interaction) poses no particular problem. In 
the symmetric phase all quantities can be computed in 
a perturbative expansion in h^. For —^ a nontriv- 
ial exact limiting solution of the functional integral has 
been found [l^. Nevertheless, for large a new strong 
interaction appears and the bosonic fluctuations become 
important. Indeed, the crossover for a broad Feshbach 
resonance amounts to a theoretical challenge. One has 
to deal with a genuinely non-perturbative setting. 

In the present paper we have included the molecule 
fluctuations by the solution for Schwinger-Dyson or gap 
equations. For T = the results are quite satisfactory 
as shown by the comparison with quantum Monte Carlo 
simulations [s^l in flg. [51 However, as T increases to- 
wards the critical temperature Tc our method becomes 
less reliable since the details of the treatment of the 
molecule fluctuations play an increasing role, cf. fig. [TT] 
The main problem results from the neglected momentum 
dependence of the four-boson- vertex A^. This is related 
to a general difficulty for Schwinger-Dyson equations: the 
momentum integrals involved require the knowledge of ef- 
fective couplings in a large momentum range. A promis- 
ing alternative may be the use of functional renormaliza- 
tion [5^ . At every renormalization step only a narrow 
momentum range plays a role and the momentum depen- 
dence of suitably chosen effective vertices is much less 



important. First results of a functional renormalization 
group treatment can be found in [ioj . 

In this light the present paper should be viewed as 
a starting point for more accurate investigations with 
further functional integral techniques. It is well suited 
for systematic extensions, among which we would like 
to stress a more appropriate treatment of the momen- 
tum dependence of the couplings, an extended inclusion 
of the effect of boson fluctuations, and the modification 
of the fermion propagator by the renormalization effects 
originating from mixed fermion-boson contributions. 

It remains to be seen if theoretical improvements, 
together with a reduction of systematic experimental 
uncertainties, will finally lead to an understanding 
of ultracold fermionic atoms with high quantitative 
precision. 
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APPENDIX A: PARTIAL BOSONIZATION AND 
PARTICLE DENSITY 

In this appendix we extend the formulation for the 
functional integral in such a way that we are able to 
discuss situations which arc (i) inhomogeneous and (ii) 
beyond the "small density approximation" (SDA). The 
treatment of inhomogeneities is particularly desirable in 
the context of ultracold gases which are prepared in traps 
of various shapes. In particular, we show how the "lo- 
cal density approximation" emerges from our formalism. 
Further we consider situations beyond small densities. 
The price to pay is the inclusion of a further functional 
integration over a now fluctuating field tr. In the SDA, 
treating the chemical potential as a source term is appro- 
priate, and we will quantify this statement here. 

We will start the discussion for a single fermion field. 
This describes a situation far off a Feshbach resonance, 
where molecules of effective collective states are unim- 
portant. Actually, the discussion in this appendix cov- 
ers a very large class of fermionic systems with effective 
pointlike interactions. Besides ultracold fermionic atoms, 
it may be applied to neutrons (e.g. in neutron stars), 
dilute gases with short range interactions (e.g. dipole 
interactions) and also covers certain aspects of electron 
gases (where the Coulomb interaction is replaced by a 
pointlike repulsion). We then extend the discussion to 
the case of fermions and bosons in order to account for 
strong interactions via a Feshbach resonance. The con- 
cept presented here will also be technically useful for the 
analysis of strongly interacting ultracold atoms in the 
frame of the functional renormalization group. 
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1. Functional integral 

For an arbitrary fermionic theory with "classical" ac- 
tion Sf [ip] we define the partition function as a functional 
of the local source J{x) 



With VFf[^^] = lnZF[J] the (relative) particle density 
becomes 



nA{x) = {'ijj^x)i^{x)) 



SWf 
5J{x) 



(A2) 



The physical particle density is related to nf^{x) by a 
constant shift that we have discussed in sect. IIIII The 
source may be composed of a bare chemical potential 
/^A and a local potential Vi [x) , 



J{x) ^ HA - Vi{x). 



(A3) 



For ultracold atoms the local potential Vi{x) represents 
the trapping potential. 

In this section we develop a general functional integral 
approach for the computation of the density n(x). For 
this purpose we introduce a field a{x) which is conju- 
gate to J{x). It will play the role of an effective chemical 
potential which differs from eq. (jA3[) unless the den- 
sity is small. Our approach can be used for arbitrary n 
and is not restricted to a small density approximation. 
The effects of density fluctuations can be incorporated 
by functional integration over a via partial bosonization. 



2. Partial bosonization 

Partial bosonization is achieved by inserting a Gaus- 
sian integral which is equivalent to one up to an irrelevant 
multiplicative constant (Hubbard-Stratonovich transfor- 
mation [13) 



Zf[J] 



VijjV&expl- SfI^^] 



{-SfI^P] 



dx 



J{x)ijj\x)ijj{x) 



-\m^U{x) - J{x) - i^^t(^)^(2,)^'l .(A4) 
2 V TO^ / J 



The partition function Zf[J] can now be computed 
from an equivalent functional integral which also involves 



The split between a constant part in Vi and fi/^ is arbitrary. For 
interacting atoms in a homogeneous situation (V; = 0), /Ia is 
related to the true chemical potential by a constant shift (see 
below). For electronic systems V; typically corresponds to an 
electrostatic potential. 



bosonic fields (t{x) 

Jdx(^J{x)a{x)-^j'^{x)^^ (A5) 

with 



(Al) SB^SF[i^]+ fdx llm^a^ - ai^^ + -^{tp^f 
J [2 2m^ 



(A6) 

Here the expectation value u{x) — {(^{x)) is related to 
the density nA{x) by 



nA{x) 



6Wf 
5J{x) 



a{x) = {a{x)) ^ —{^\x)^P{x)) + J(x). (A7) 

After partial bosonization the new action 5*^ contains 
a mass term for ct, a Yukawa interaction aip^il^ and 
an additional four-fermion interaction ^ {ip^ ip)'^ / {2m'^) . 
The resulting explicit four-fcrmion vertex in Sb therefore 
becomes 



Av, — A + 



(A8) 



We can choose such that A^ vanishes by requiring 



(A9) 



The limit of non-interacting fermions obtains then for 
— + oo. We note that such a cancellation is possible 
only for A < 0. It is not compulsory for our formalism 
an we will not always assume the relation (jA9p in the 
following. 

The partition function Zf (|A5|) . (j A6[) is closely related 
to the standard formulation of a scalar-fermion-model 
with 

Zb^ J Vi;Vaexp^^-SB[^,ij]+ Jdxj{x)a{x)^ (AlO) 
by {WB=lnZB) 



j{x) = m^J{x), Wb^Wf + — JdxJ^ix). (All) 

In this formulation we may define the (one particle irre- 
ducible) effective action F by the usual Legendre trans- 
formation 



Wb + dxj{x)a{x), a{x) 



6Wb 
Sj{x)' 



(A12) 



For a given source j{x) the expectation value cr{x) obeys 
the field equation 



_SI_ 

Sa{x) 



(A13) 
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The effective action may be decomposed into a classical 
part r^; [cr] = Sb [o" = cr, tp = 0] and a quantum part Fg 



by the physical particle number density nx + n. This 
yields 



T = r,i+Tg='-^ Jdxa^{x)+rg[a]. (A14) N ^ J d^xiuA+n) ^ fri^exp (^^^^^^ +1^ \(A22) 



We next turn to the relative particle density n\{x). It 
can be computed from F by decomposing the field equa- 
tion 



a{x) = J[x) + <7i{x), 



5T 
5a{x) 



as 



m'^a{x) 



nA{x)=m''ai{x) = -^^(^J{x)+ai{x)y (A16) 

This establishes the exact general relation between nA{x) 
and the a - functional derivative of the quantum part of 
the effective action. 

The limit |cti| ^ |J| corresponds (for fixed J{x)) to 
small densities or small interactions {m'^ oo) according 
to 



da{x) 



= m^J{x) (A15) 



\nA{x)\ < m'^\J{x)\. 



(A17) 



In this limit we may expand 



Sa{x) 
with 



^(jix)+ai{x)) = -no{x) + b{x)ai{x) + ... (A18) 
X) \ / 



no{x) 



This yields 



Sa{x) 



J{x) 



b{x) 



S^a{x) 



J(x) 



njy{x) 



no{x) 



1 + h{x)/m? 



(A19) 



(A20) 



We emphasize, however, that eq. (|A16p remains valid for 
arbitrary densities. 

The explicit computation of {x) needs an evaluation 
of Tq. In mean field theory the fluctuation part F^ is esti- 
mated by including only the fermionic fluctuations in the 
functional integral (jAlOp . while keeping a{x) = (t{x) as a 
fixed "background" . In this scheme, the generalization of 
eq. pTjl for the relative particle density reads, according 
to eq. (|XT6)) 



d^xn^ix) = — 



Sa{x) 



-— Tr tanh 
2 



2T 
(A21) 

where the remaining trace Tr is over three dimensional 
phase space and spinor indices a - in momentum space 
it reads Tr=V3 / ^jjryr J2a^ v/ith V3 the volume of (three- 
dimensional) space. At this point the construction from 
sect, nira can be performed, replacing the relative (n\) 



This formula has a simple interpretation. The trace Tr 
over an operator A can be evaluated in a complete or- 
thonormal system of complex functions /,„ (x) , 



^ *^/rn' (*^) (*^) — ^ra'mi 

d^xf„-^,{x)Afjyi{x) — A^i^ 

Try4 = ^ ^ Aram ■ 



(A23) 



With a{x) = n — V{x) (in analogy to eq. (|A3j)) we can 
define the Hamilton operator H — Pp + V (with P — 
— A/(2M) for nonrelativistic atoms). Choosing the fm 
to be eigenfunctions of the Hamiltonian with eigenvalue 
Em eq. (|A22p becomes 



N = 



E 



exp 



Em — i» 

T 



1 



(A24) 



and we recognize the well known fermionic occupation 
number. In the limit of noninteracting atoms (m^ — *■ 
00) we can cut the Taylor expansion in eq. (jA18|) at 
the lowest term, h\ = uq. In consequence, STq/6a is 
evaluated with the local potential Vi in eq. (|A3p . In 
this limit one has a — fi ^ Vi and V therefore equals the 
trap potential V/. As it should be the Hamiltonian H 
reduces to the quantum Hamiltonian of a single atom in 
a potential and the density becomes independent of m^. 



3. Universal field equation for the density 

The crucial advantage of our formalism is the possibil- 
ity to compute the effective action T[a] without any spec- 
ification of the local potential (trapping potential) Vi{x). 
The detailed geometry of the trap only enters at a second 
step when one solves the field equation (jAlSp . This offers 
the great advantage that the fluctuation problem can be 
solved in a homogeneous setting, i.e. one encounters stan- 
dard momentum integrals in the loop expressions rather 
than summations over spherical harmonics or other func- 
tion systems adapted to the geometry of the trap. The 
results are universal and apply to all geometries, provided 
the inhomogeneity is sufhciently weak. 

For <7{x) mildly varying in space and independent of r 
one may use a derivative expansion 



r[a]^ ldx'>Uai(7) + -A,ia)VaVa 



(A25) 



We emphasize that the derivative expansion always ap- 
plies for situations close enough to homogeneity. The 
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following results are therefore valid independently of the 
way how T[a] is computed - in particular, they do not 
rely on MFT. However, for definiteness, let us give the 
result for in the latter scheme {U{a) can be read off 
from eq. (jl9p for 7^ = 7). For this purpose, we proceed 
along the lines in sect. IIII CI and extract from Tq,MFT 
the term quadratic in Scjq 



(A26) 



PAQ) = I { 



1 



{Ppm - 'y){PFiQ' + Q) - <7) 

1 



{PpiQ') ~ a)iPF{Q' ^ Q) 



with Jq, = X^n'^/ (2it)^ • Inserting the background 
(t{X) = (t + cxp{iQX)SaQ + c:xjp{—iQX)6aQ into eq. 
(|M5l) yields 



A„{a) = lim 



d 



(A27) 



1 f cPq ( tanh 7 2 cosh^ 7 — 8 

JmT^J (27r)3 \ cosh^ 7 " 9Afr coshS 



Here we have chosen the definition of U such that the 
physical particle density from eq. (|23p is reproduced in 
eq. (lAaOl) . 

In the second line of eq. (|A30[) appears a shifted or 
"additively renormalized chemical potential" , 



Ai = MA 2- 



(A32) 



This combination is the true "physical" chemical po- 
tential as can be seen along the lines of eqs. (|22|23p . 
Indeed, the interaction term in the operator language 
cx (a^(x)a(a;))^ is translated to the interaction term in 
the functional integral (|A1|) for the microscopic action 
by the use of 



2m^ 2m'' \ 



(A33) 



The term {h/m,'^)'ip^{x)ip{x) precisely shifts /i to /2a in 
the functional integral (jAl[) . A detailed account for a 
similar "chemical potential shift" for bosonic atoms can 
be found in f6^. 

Hence the second part of eq. (jA30|l yields a linear 
relation between the physical particle density n and a, 



a = ^i-Vi^ ^. 



(A34) 



We find that Acr{(j) is ultraviolet finite. The zero tem- 
perature limit of this expression reads (for cr > 0) 



A„ = 



1 2M 



127r2 V cr 



(A28) 



and diverges for ct ^ 0. 



a. Density field equation 

We are interested in time-independent situations and 
therefore consider a{x) = u{x) independent of r. Varia- 
tion with respect to a yields the stationary field equation 
(with U'„ ^ dUa/da etc.) 



UU'T)-A4a)Aa--A'Aa)Vc7Va = m^flA~Vi) (A29) 

where a and Vi depend on x. Once a{x) is found one can 
compute the particle density from 

1.' 



m (a - ^A + Vi -\ 



U — Ua — -rri^a^ — ha. 



(A30) 
(A31) 



The prefactor has been determined numerically with an error of 
less than 0.3%. 



This will be our central equation relating the density n, 
the effective chemical potential a, the chemical potential 
fi and the local trap potential VJ. We observe that a 
great part of the "ultraviolet divergencies" (for A — > 00) 
present in the functional integral description is related 
to p,\ and h. These divergencies do not appear if the 
"physical quantities" fi and n{x) are used. 

Eq. (|A34[) can now be used to eliminate a in favor of 
the physical particle density. Insertion into the first eq. 
(|A29p yields the central field equation for the density in 
an inhomogeneous situation, 

n - -^A„ (^-Vi + At) An 



Vn Vn 



+ ^A'A^,-Vl + A,)[^V^ 



T) 11 

-U'ifi -Vi + ^)- A^pi -Vi + -^)AVi 



1 ft 

+ -A'^{^ji-Vi^^)'^Vi'^Vi. 



(A35) 



Once the functions U'{a) and A^{(j) have been computed 
this equation permits the determination of n{x) as a func- 
tion of /i. For its solution the boundary conditions of n{x) 
have to be specified appropriately - for trapped atoms 
n{x) has to vanish away from the trap. In practice, one 
often does not know /i in a given experimental situation, 
but rather has information about the total particle num- 
ber N — J (fxn{x). In this case one may compute N{^) 
by integrating the solution of eq. (|A35|) . Our setting 
can be reformulated in a more intuitive form in terms of 
a "density functional" r[n] to which we proceed in the 
next section. 



30 



b. Density Functional 

We can write the effective action in a more intuitive 
form as a functional of the particle density. It is con- 
venient to add the source explicitly such that the field 
equation becomes 



5T ^ 

Sn{x) 

Using eq. (|A34p one has 

r[n] = r[cr] + Jdxlm^iVi - n) - n](T 

t2 



(A36) 



(A37) 



= jdx\u{ 



+ ^i-Vl) + 



n 
2m? 
n 



where again a density independent term has been 
dropped. The explicit form of eq. (jA36[) is given by 
eq. (fM5l) . 

In particular, a homogeneous situation is governed by 
the effective potential 



2tv? 



U{n). 



(A38) 



In this case we can use n as the independent variable, 
replacing ^. In turn, the chemical potential follows 
from the minimum condition dUn/dn — 0. In practice, 
one first solves for a{n) by inverting 



U'ia) 



(A39) 



In MFT this step does not depend on m?. The inter- 
action strength only enters the determination of /i 
through eq. (jA34[) . Similarly, in the local density approx- 
imation we can trade /x for the density at a given refer- 
ence location xq, uq = n{xo). Using again eq. (|A39[) . the 
computation of the density profile employs 

a{x)^a{no) + Vi{xo)~Vi{x) + '^^^^^. (A40) 

For a given central density no (e.g. xo — 0) and given 
trap potential Vi (x) — Vi (0) we can now evaluate U' as a 
function of n and determine n{x) from U'{a{n)) = —n. 

On the other hand, if n{x) is known we may use a 
functional of the variable a for fixed n[x). Using 



5T_ 



m a — n + n 



(A41) 



we may define 

f = r + y dx^l^n — nja 2"''^^}' 



sr 



(A42) 



Here n{x) is considered as a fixed function. The corre- 
sponding potential 



U - —ct'^ + {n-n)a = U + n(T (A43) 

is particularly useful for a homogeneous situation where 
is a fixed constant. The solution for cr(n) corresponds 
to the minimum of U . 



c. Small and local density approximations 



by 



In the small density approximation (SDA) a is given 



cr = ^ - Vi{x). 



(A44) 



We can now specify the validity of the small density limit 
more precisely and replace the condition (jA17P by 



n < m^\n - Vi\. 



(A45) 



In lowest order in the small density approximation one 
obtains the density of non-interacting fermions in a local 
potential as 



-u'(ji-Vi)- A„{^l-Vl)^Vl 
-]^A'^{^,-Vl)^Vl^Vl. 



(A46) 



We emphasize that eq. (|A35I) and its low density limit 
(IA46|) are rather universal formulae which determine the 
density of fermions in a local potential. Their possible 
applications reach far beyond the particular problem of 
ultracold trapped atoms. 

Another useful approximation, the local density ap- 
proximation (LDA), obtains by setting = in eq. 

n = -[/'(^; + _), (A47) 
where we define the "local chemical potential" 



fii{x) = /i - Vi{x). 



(A48) 



Its validity does not require a small density but rather 
that the derivative terms in eq. (|A35[) (or (jA25|) ) are 
small as compared to U'. For a given size of A„ this 
always applies for a sufficiently homogeneous trap. In- 
deed, the local character of eq. (|A47p guarantees that 
weak changes in /i;(x) result in weak changes of n{x). 
The error of the LDA can now easily be estimated by an 
iterative procedure. Inserting the solution of eq. (|A47[) 
into eq. (jA35P allows for a simple direct computation of 
the subsequent terms. Obviously, the error depends on 
the size of Aa- which is often a rather small quantity (see 
below) . 

If both the density is sufficiently small and the trap 
sufficiently homogeneous one can work with the small 
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local density approximation (SLDA). This results in the 
simple formulae 



dU_ 
9^ 



and 



a = ^ii. 



(A49) 



(A50) 



In this approximation we can regard cr(a;) as a fixed ex- 
ternal parameter and compute n{x) by eq. (|A49|) . We 
will find that for realistic ultracold atom systems like ^Li 
the LDA is valid whereas eqs. (IA49[) . (jA50|l are oversim- 
plifications (except for the BEC regime far away from the 
Feshbach resonance). 

We emphasize again that our method is valid quite 
generally and not bound to the case of small density. If 
necessary, mean field theory can be improved by includ- 
ing the bosonic fluctuations in the computation of U and 
Aa- by performing the functional integral over a. Our 
method can also be applied in the presence of additional 
condensate fields. As a simple application we compute 
in appendix |B] the density at low T in the mean field 
approximation. In particular, this shows that the dilute 
gas of ultracold atoms is a metastable state, the ground 
state being a liquid or solid. 



have to be fixed by appropriate observable "renormal- 
ized" parameters. We have also included a possible local 
self- interaction of the molecules ~ and between 

free atoms and molecules ~ 2(3 /w?'. We concentrate on 
Vm = 2V; and (i = $ = 1. This simply counts the 
molecules as two atoms as far as the local interactions are 
concerned, e.g. the local self-interaction is ~ {nF+^nMY' 
and the energy in the trap potential is rip + 2nM- We 
note that in this particular case h drops out. (The cor- 
rections to /i are computed similar to eqs. (jA32|A33|) .) 

Replacing in the Hubbard-Stratonovich transforma- 
tion (jA4p tp'^ip -0'!'^ -t- 2(/)*0 all steps in sec. lA 31 can be 
performed in complete analogy, with the only difference 
that Zb in (jA10[) involves now also an integration over 
(j> and an appropriate source for (f>, and Sb reads (for 



Se 



dx 



A 
4M 



i^A - 2o-) 



(A54) 



This is a simple model for fcrmions with Yukawa coupling 
to scalar fields 6 and a. 



Fermions and Bosons 



We can proceed in complete analogy to sects. lA 1[ lA 21 
with Zp replaced by 

Zfm[J] = JvipV^e^pl^-SpASJ] (A51) 
Jix){i^^{x)iP{x) + 2$*ix)$ix))]}, 



+ dx 



and 
S In Zp M 



n{x), J{x) ^ n+i^-P) — -Viix). iA52) 



5J{x) 

The classical action reads now 



SpM = I dx 



+'^* (^-- ^ + + 5 (1 - + /3) + Vm {x) - 2Vi (x) ) $ 



2(3 



4M m 

2^2 2/3 



(V>t^)(^*<?!)- -^{<|)*^P^e'l|J-<^'l|Jh^p*) 



(A53) 



and the partition function describes the coupled sys- 
tem of atoms in a local potential V/ and molecules in 
a "molecule potential" Vm- (The term cx — 2V/0*(?!) in 
SpM cancels the corresponding term from —2J(j)*(f).) The 
atoms are coupled to the molecules by the Yukawa cou- 
pling h^. Again, the bare parameters j^a, and 



a. Effective action 

The effective action T[a, cj)] obtains again by a Legendre 
transform similar to eq. (jA12[) and obeys now 



2 2 

— = m til = m (J — n. 

0(7 



(A55) 



In particular, a homogeneous situation is characterized 
by an extremum oiU = Ua — (m^/2)(T^ + ntr. 



Ana 



■<l)*4> - 2a(j)*(j) + Ui^U + na. (A56) 



Similar to sect. (|VD[) we may proceed beyond the Gaus- 
sian functional integration for ip (MFT) by adding to 
Ui a piece from the one loop contribution from the 0- 
fluctuations 



(A57) 



The contributions from dressed uncondensed molecules 
is now given by 



da 



= -2n 



M 



(A58) 



Indeed, eq. (jASSp follows from eq. ([51)1 if /i is replaced in 
the classical bosonic propagator by the effective chemical 
potential cr. 
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6. Field equations 

Collecting the different pieces the extremum of U (eq. 
(|A56|) 1 occurs for 



-r— = = ?i — 20 (h — npo — 2nM- 
da 



(A59) 



This is precisely the relation (|115p as it should be. (Recall 
np.o — —du[^^ /da.) In other words, the density obeys 



dU_ 



(A60) 



With J — fj, we may derive the relation between the 
effective chemical potential a and the chemical potential 
A* 



II = a — n/rn^ 
directly from the identity 



(A61) 



J ViPVaV4>-^ exp {- SB + rn^ j Ja] ^ 0. 

(A62) 

The evaluation of ?7(ct, (j)) is now sufficient for the compu- 
tation of the total atom density n. For the homogeneous 
setting one can therefore determine a by 



dii 1 
da 37r2 

c. Effective chemical potential 



(A63) 



In summary, our problem is now formulated as a func- 
tional integral for a Yukawa model. In the small density 
approximation we can treat a{x) = fi + Vi{x) as an ex- 
ternal parameter. The partition function becomes 

Z = y Vipvlj) exp [ — 5*3 [^A, 0; cr] -I- source terms](A64) 

where Sb is given by eq. (|A54[) and the density obtains 
from eq. (IA63|) . Beyond the small density approximation 
a is treated as a field and the partition function involves 
an additional integration over a. In the limit where the a 
- fluctuations can be neglected we may consider the effec- 
tive chemical potential a instead of /i as a free "external 
parameter" . In particular, this offers for the homoge- 
neous case the advantage that we are not bound to the 
validity of the small density approximation (which may 
not be accurate for realistic systems as ^Li). It is suffi- 
cient to compute the value of a{T) for a given density or 
given kp. 

For the homogeneous case we can both deal with situ- 
ations at fixed effective chemical potential a or at fixed 
density n. For a fixed a we may choose an arbitrary fidu- 
cial kp (not determined by n) and do all the rescalings 



with kp (instead olkp). One can then work with a fixed 
value of a and compute n or kp by the relation 



du 



1 



1.3 
l^p 



37r2 fc3 • 



(A65) 



Many experimental settings can be idealized, however, 
by a fixed value of n. Then the choice kp = kp and the 
determination of a via eq. (|A63p provides directly all 
results for fixed n. For an inhomogeneous setting it is 
sensible to use a suitable fiducial kp. 



APPENDIX B: METASTABILITY 

At low T the thermal equilibrium state of atoms is a 
liquid or solid, with density essentially determined by 
the "size" of the atoms (typically set by the van der 
Waals length). When we deal with a dilute gas of atoms 
at ultracold temperature we obviously do not consider 
the stable thermal equilibrium state which minimizes the 
free energy. Indeed, we will see in this section that the 
metastability can be captured within our formalism if the 
weak cutoff dependence 0{A~^) is not neglected. For this 
purpose, we consider a homogeneous system and neglect 
the effects from molecules. 

Under the above circumstances the field equation 
(|A35P reduces to a simple self-consistency relation for 
the density, 

n = -U' f/i + 

d^q 1 



= 2 



(27r)3 e('3V(2M)-^-n/m2)/T ^ I 



(Bl) 



For the gross features we first consider the T — > limit 
where the last equation reduces to 

The resulting cubic equation for y — n/dim^) 



(B3) 



may have several solutions, depending on k. We concen- 
trate on /i > and consider first k ^ 1. The solution 
with small y ~ k^^/^. 



37r' 



:(2M^) 



3/2 



(B4) 



describes a dilute gas of atoms. In lowest order the re- 
lation between n and ji is indeed independent of rn? and 
A. For large n the second solution with positive y occurs 

for y fa K 



8M3 



n2 



(B5) 
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while the third solution has negative y and should be dis- 
carded. As K decreases (i.e. /i increases) the two solutions 
approach each other and melt for Kc — 27 / 4 , yc — 2 or 



2/iTO 



3M3 ' 



(B6) 



which delimits the metastable gas phase. No solution 
with y > exists for k < Kc- The stability of the solution 
depends on the second derivative of C/„ with respect to 
n, i.e. for T ~ 



1 



1 - 



112- 



(B7) 



APPENDIX C: BOSON PROPAGATOR IN MFT 

In mean field theory the gradient terms for the fields 
a and (/> obtain by evaluating the fermion loop in an 
inhomogeneous background field. We write the inverse 
fermion propagator as V+T with V the part independent 
of a and (j). We decompose the fields in a position inde- 
pendent part and a harmonic (single momentum mode) 
small local perturbation, 



T{X, Y) - Th + T^nh 




The solution rii (with small y) turns out to be stable 
{d'^Un/ dv? > 0) whereas the solution n2 (with y « k) is 
unstable. For Kc one has d^Un/dn? = 0. 

What happens for n > ric^ In order to study this ques- 
tion we cannot neglect the ultraviolet cutoff A anymore. 
In the evaluation of eq. (jBip the upper bound of the 
momentum space integral is actually given by qmax = A 
instead of infinity. For A < {2M{fi + n/m^))^/^ this 
multiplies the r.h.s of eq. (|B2p by an additional factor 
A^/{2M{fi + n/TO^))^/^ such that 



A^ n 
37r^ TO^ 



A2 

^i>—■ (B8) 
^ 2M ^ ' 



This solution is again stable with d^Un/dn^ = l/m^. We 
associate it with the liquid or solid thermal equilibrium 
state. Indeed, the free energy for Ugs = A'^/(37r^) is much 
lower than for the solution ni . Also the density is deter- 
mined by the effective size of the atoms ~ A^^. Inversely, 
we can define our cutoff by the density ngs of the liquid 
or sohd ground state for T = 0, A = (37r^ngs)"^^^- Of 
course, in this region our approximations are no longer 
valid, but the detailed properties of the liquid or solid 
state are not the purpose of this paper. 

Beyond the limit T — > 0, a computation of the phase 
boundary for the existence of a metastable gas for arbi- 
trary temperature requires the solution of the condition 
for the melting of the stable and unstable extrema, 



^ 



= — 



with 



m 



U" 

— / — ^ cosh~^ 7. 
2T / 27r3 ' 



(B9) 



(BIO) 



One of the dependent variables, /i or n, has to be elim- 
inated via the equation of state I|B1[) . For the range of 
T considered in this paper the influence of the temper- 
ature is negligible and the qualitative result (jBSp is not 
affected. 



where ea/3 is the 2 dimensional totally antisymmetric 
symbol and 

4>{X) = ^ + S4>exp{iKX), (C2) 
a{X) = <T + Sa exp{iKX) + Sa* exp{-iKX). 

The one-loop fermionic fluctuations can now be expanded 
around the homogeneous fields. 



r(MFT)^_lTnn(p + J-,+J-,„) 



(C3) 



^-^TT\ii{V+Th) - ^Trln(l+ T,a{V+Th)-') 
= VUmft + \Mmn{r + + 0{Tt^). 

The terms with odd powers in Tinh vanish due to trans- 
lation invariance. 

For a practical computation, we switch to Fourier space 
where 



(C4) 



and 



7^(Ql 



{V + Th)-\Qi.Q2) = Tl{Qi)5{Qi - Q2), (C5) 

\-5af3{PF{Ql) - ^) -^a(3h^(t>* ) 

1 

[Pf{Qi) - a][PF{~Qi) - a] + hl4>*^' 

We concentrate first on the gradient term for the 
molecules {A^) and set 5a — Q. Then the inhomogeneous 
part reads in momentum space 



^in{Ql, Q2) — hcl,eap 



*SQ1:Q2+K 



-K 



(C6) 
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{6r.s = S{R - S)). This yields 
1 



^Tii[T.,niV + Th)-^]^) (C7) 
\tT J T,„iQi,Q2)n{Q2)T,niQ2,Qi)niQi) 



Qi 



where tr is the trace over the "internal" 4x4 matrix. 
One obtains 



^MFT _ 



(C8) 



PFhQ')-<T 



[PpiQ') ~ a][PFi-Q') - a] + hlq 
PF{Q'-K)-a 



[Pf{K - Q') - a][PF{Q' -K)-a]+ /i^^*^' 
Insertion of the ansatz (IC2I1 into the effective action 



(C9) 



yields eq. p6|) . We note the simpler form for = 
(symmetric phase) where 



V^{K) = 27rinT + 



g 

4M 
hi 



(CIO) 



Q 



[PF{Q)-a][PF{K-Q)-a\ 



The quantum corrections differ from P^ (eq. (|A26[) ) by 
the overall factor — ft,^/2 and the different momentum 
structure in the denominator. 

We note that in the superfluid phase the symmetries 
would also be consistent with a gradient term of the form 
/g (j)*5(j){Q)(j)*5(j)(~Q). The coefficient of this term can- 
not be computed with the ansatz (|C2[) - it would re- 
quire a more general ansatz 4'{X) — (l)-\-d(j)+ exp{iK X) + 
6(j)- exp{—iKX). We will omit this term in the present 
paper. 



1. Momentum dependence 

For spacelike momenta the loop integral depends on 
the square of the external momentum only and we define 
Aci,{q) and by eqs. (|41l42p . The fluctuation correction 



Note the slight abuse of conventions. Here, V stands for the 11 - 
entry of the inverse propagator matrix in the </>*, {j> basis, instead 
of the full matrix. 



to at n = is given by 



/(?') 



3 [((2m-f l)7rT)2 + /2(g')+r] 
/(?' - q) 



[{{2m+l)7rT)^ + P{q' -q)+r] 
~K f d^q' f{q')f{q' - q) 



(CU) 



2 J (2^)3 f{q') - f{q' - q) 

r 1 , Vb 1 , y^) 

X < —= tanh — = tanh — — > 

IVb 2T 2Ti 



with f{q) = qy2M-a, a = f{q')+r, b = p{q'-q) + r. 
The momentum dependent gradient coefficient reads 



Mq 



2, _ 1 AP^jq^) - AP^jO) 
AM q^ 

- Af^^AA,{q\ 



(C12) 



As argued in sects. |V] and IVIl the physically relevant 
quantity is A^/Z^. We have plotted the momentum 
dependence of A^{q)/Z^ in the broad resonance limit 
/i^ ^ 1 in fig. 01 At large momenta, the gradient co- 
efficient slowly tends to zero. This has no impact on 
observables, since the thermal distribution functions are 
suppressed at lower momenta already (see below, also fig. 
[SJ . It might be an artefact of the neglected momentum 
dependence of Z^. 

In the present paper we neglect for the numerical re- 
sults the momentum dependence of A^ and approximate 
A^{q) = Acf, = A^{q — 0). Let us discuss here the valid- 
ity of this approximation in the broad resonance regime 
at the critical temperature. The impact of the momen- 
tum dependence of A^{q)/Z^ is most clearly seen for the 
Bose distribution in fig. [5l There we explicitly com- 
pare the results with a momentum dependent ^^(q)/^^ 
and an approximation of constant A^/Z^. In the BEC 
regime the error is very small with A^{q)/Z^ close to the 
classical value 1/2 for all q. The main difference from 
the classical propagator concerns the renormalization of 
the mass term m^. In the crossover regime the approx- 
imation of constant A^/Z^ underestimates the number 
of molecules with large q^ , but the error remains accept- 
able. In contrast, the deviation from the result with the 
classical molecule propagator is already substantial. In 
the BCS regime the underestimate of NM{q) for the dom- 
inant range in q is quite substantial. Though the over- 
all effect of the boson fluctuations is small in the BEC 
regime, this may affect the quantitative results for the 
number density of molecules and the condensate. In view 
of flg. [5] (c) the estimates of CIm and flc and f^M, in 
the present paper are most likely too small. 

We next discuss A^ = A^{q = 0) more explicitly. With 
A^ = 1 /AM + AA^ we can compute AA^ as the term 
linear in q^ in the Taylor expansion of AP^ (|C11[) . Using 
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A, 



2MArh we find the result 



2TrinTf(q'). This yields {uj = 2TmT) 



17 



A, 



2 ' 288f3y (27r)3 
X [tanh 70 - 70 cosh~^ 70] + 27^70(7^ 
X [70 cosh"'* 70-6 tanh 70 - 270 tanh^ 70] , (C13) 



simplifying in the symmetric phase 70 = 7 to 



2 48f3y (27r) 



tanh 7 



7 cosh 7] . 
(C14) 



The loop correction to the gradient coefficient is strictly 
positive and monotonically growing for a increasing from 
negative values (BEC side) to its saturation value a = 1 
on the BCS side. For the BEC regime it vanishes in 
the limit a — > —00. However, the physical quantity 
Aif, =A^/Ztp approaches the finite value 1/2 in the BEC 
regime - this is an indicator of the emergence of an effec- 
tive bosonic theory. For a ~ 1 instead, A^jZ^ is much 
larger than the value for elementary pointlike bosons, 
1/2. Indeed the integral is dominated by modes peaked 
around a, explaining the strong increase as 17 — > 1. 

In the limiting BEC and BCS cases, A^p can be ap- 
proximated by much simpler formulae. The BEC result is 
given in app. [D] eq. (|D3[) . In the BCS regime, the critical 
temperature is very small (T < 0.1) and 0.4 ct < 1. We 
then find an approximate behavior for the f-dependent 
gradient coefhcient 



1 7C(3) hld^'^ 

2 127r4 4^2 -f f 



(C15) 



In the temperature range of interest the classical contri- 
bution to A^ is small (^0^''' = 1/2). Neglecting it and 
restricting to the symmetric phase {(j) = 0}, this is con- 
sistent with the symmetric BCS result in 11[, and we see 
how the condensate regulates the divergence for T ^ 0. 



ImAp0(a;,g = G) = c^/j^T^ j ^f{q') (C16) 

x{((2m + l)7rT)2 + /2(g')+r}-i 
x{((2m + 1)ttT + ujf + f{q') + r}-\ 



Obviously ImAP0 vanishes for n = {oj = 0). For n 7^ 
we can perform the Matsubara sum (we suppress the 
argument of /), 



UJhj f <Pq' tanh(v//2 -f r/2r) 
(27r)3 (c^/2)2 + /2+r ' 



(C17) 



We may define the coefficient of the Matsubara frequen- 
cies as (wi = 2-kT) 



^0,r(wi) 



Im(^0(a;,O)) 



Im 



LO 

^0(^1,0) 



(CIS) 



K f d^'q' tanh(/F+^/2r) 
47 (27r)3 {-KTY + p+r ' 



We can study the bosonic propagator as a function 
of m. The loop corrected imaginary part of the inverse 
boson propagator can be brought into the form 

ia;,„(l -f c(g^, T, cr, m^)), ujm = 2Trm. (C19) 

Each Matsubara mode is renormalized by an m- 
dependent quantity. In the present paper we neglect 
these corrections, i.e. we take c{q^,T,a,m'^) = 0. 



APPENDIX D: ANALYTICAL RESULTS IN THE 
BEC REGIME 



2. Frequency dependence 

Similarly the loop correction can be evaluated as a 
function of external Matsubara frequency ujm {Q = 
(Lu„i,q)) for vanishing external spacelike momentum q. 
This amounts to the renormalization correction to the 
operator (j)*dr<j)- In momentum space this corresponds to 
the part ~ \0Jn in the inverse molecule propagator. Hence 
we only need to consider the imaginary part of the loop 
integral. The denominator in the integrand in eq. (jCSp 
is real and the imaginary part of the numerator becomes 



We exploit the fact that in the BEC regime a/2T 
—00, which means that we can replace the functions 
tanh 7 l,7C0sh~^7 ^ 0. In the superfluid phase, 
we use additionally r/|tT| — > 0. The loop integrals can 
then be evaluated analytically. They are temperature 
independent. Furthermore, their values coincide in the 
symmetric and superfluid phase, if terms of 0{f /a) or 



As in eq. I|C9| | , V stands for the 11 - entry of the inverse propa- 
gator matrix in the (j>* , basis 
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higher order in ct are ignored. We find 



-a 



8^ ' 

H 

12877 ' 



(Dl) 
(D2) 

(D3) 
(D4) 



For the fermionic particle density contribution we find a 
term 0{f / —a), 



(D5) 



The BCS gap equation is solved by using (|D1 



c"^ = V^, (D6) 

independently of the value of f. Hence in the BEC limit, 
the relation between binding energy and scattering length 
[m is independent of the density scale kp- Indeed, many 
body effects should be unimportant in this regime. Ap- 
proaching the resonance, the impact of the pairing gap r 
becomes important and a ^ sm/^f- 

In the limit of large Yukawa couplings, we can then 
evaluate the gradient coefhcient of the effective Bose dis- 
tribution (cf. eqs. pi3lll3p V 



A, 



A^ 

Zrh 



(D7) 



Similarly, the fermionic contribution to the four-boson 
coupling evaluates to 



72 

^0 



(D8) 

Using the relation Xp — Airap/Mp between coupling 
strength and scattering length, the molecular scattering 
length in the BEC limit is given by 



(D9) 



where we have changed back to dimensionful quantities. 
This reproduces the Born approximation for the scatter- 
ing of particles of mass Mp = 2M. In approaching the 
resonance for the fermionic scattering length (crossover 
0, the bosonic scattering length however 



regime; c 



remains finite. Note that both A^ and A^'^'' are effectively 
independent of in the broad resonance limit. 



APPENDIX E: SCHWINGER-DYSON 
EQUATIONS FOR THE MOLECULE 
COUPLINGS 

In this appendix we provide details of our computation 
of the effective molecule- molecule interaction A,*. We 



work in dimensionless renormalized units. In the sym- 
metric phase we expand the effective potential u and the 
mean field effective potential umft in powers of p, 



u = mlp+^X^p'^ + ... 

UMFT - P+-^X^ P +■■■ 

The contribution from the boson loop reads 



(El) 



Ml - m, p+ -A, p + .. 



1 

AF)2 



(E2) 



= K-m^-^>+^(A,-Anp^ + 



(see below) . In the symmetric phase we evaluate u\ in 
an approximation where we truncate in quadratic order 
in p (|E2l) . 

2 (-F)2 , (B)2 



A. 



92 [/ 



9p2 



p=0 



Af^+Ai^). 



(E3) 



We determine the coupling A^ from the "Schwinger- 
Dyson" equation 



Aw) — A 



p=0 



dp^ 

3Af^A^ fd^ 



(E4) 



+2a sinh~2 a] 
a(^)+A^-/a. 



"^[(exp2a-l) ^ 



which has the solution 



Xd, — 



1-/a 



(E5) 



For m| the last term in eq. (|E4[) becomes infrared 
divergent. Divergences of this type of quantum correc- 
tions to quartic couplings are familiar from quantum field 
theory and statistical physics of critical phenomena. In- 
deed, the point = corresponds to the critical line 
(or hypersurface) for the phase transition to superfluidity 
- for negative the symmetric phase becomes unstable. 
The remedy to this infrared problem has been well un- 
derstood from the solution of functional renormalization 
group equations: the stroiig fluctuation effects drive A^ 
to zero at the critical line [s^, [H, [6l| . 

Our gap equations recover this important feature in 
a direct way. As approaches zero the negative last 
term in eq. (jE4p becomes more and more important as 
compared to A^ on the left hand side. The solution to 
eq. (|E4p impHes 



lim X^{m^) 



0. 



(E6) 
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For small values of m| in the vicinity of the phase tran- 
sition we can expand the integral in eq. (|E4[) as 

/. ^ -15f ) J + miy'- (E7) 

One infers cx according to 

h = ^A'JW (E8) 

In the superfluid phase, we expand the effective poten- 
tial around po 

u = ml{p-po) + ^{p~p,f. (E9) 

We choose again a basis of real renormalized fields , 02 
according to 

-^((/.i-Ki02), P=l{cf>l + cl>l). (ElO) 

Without loss of generality we may consider a background 
of real 0, i.e. 0i □ = 2po, (l)2,o = 0. With = (f)i - (/!)i,o 
the potential (|E9p becomes 

" = ^™0'/'2 + ^{"^1 + 2A0po)0f + \4,\/po/2(j)[(l)l 

(Ell) 

We can associate and A^ with the terms quadratic 
and quartic in 02. The dots denote cubic and quintic 
terms oc ^'f^, 020ii '/'2'^i^ that will not contribute in our 
approximation and we neglect terms of O(0^). We use 
the Schwinger-Dyson equations for the and (\)\ vertices 
which result in eq. (|143p and 

> AF) 'i^P^4>id^q_ _3r/ n2. ^-l 

A0=A^ (exp2a^-l) 

+2(a + K/2)^a0sinh~^a0]. (E12) 

Again, we observe in eq. (|E12|1 the appearance of in- 
frared divergences in the contribution oc A^ A^ from the 
Goldstone fluctuations. If we would define \^{p) by the 
02 - vertex evaluated at some value p > pq they would 
be regulated. Taking the limit p ^ p^ we obtain similar 
to eq. (|E6p \^{p po) 0. The Goldstone boson fluc- 
tuations renormalize A^ to zero, as found in (stI. IssI. [61|| . 
In our approximation A^ vanishes for all T < Tc in the 
superfluid phase. In consequence, the mass term 2X^po 
of the radial mode vanishes for all T < Tc- However, 
as we have discussed in the main text, this vanishing of 
A^ concerns only the effective vertex at zero external mo- 
mentum, whereas the loop integrals often involve vertices 
at nonzero momentum. 

Furthermore, the contribution from the fluctuations of 
the radial mode in eqs. (|E4p and (|E12p are not treated 



very accurately. First, the 0']^0| vertex contains in prin- 
ciple a contribution oc f^po {i^cf coefficient of the con- 
tribution oc {p~ pqY) which shifts A^^^ + 2iy''p po 
and is neglected here. Second, the structure of the in- 
verse propagator of the radial mode is actually not sim- 
ply A^q'^ with constant A^. Indeed, the effective quartic 
coupling A^ only vanishes for zero external momentum, 
with a typical momentum dependence A^ oc \q\. For a 
definition of a mass term at q = this is consistent with 
A^po 0. However, this effect will then become visi- 
ble as an infrared divergence of the gradient coefficient 
for the radial mode, A^,r oc pol?!^^ (which differs from 
A^ for the Goldstone mode). In consequence, one has 
P,p{q 0) oc \q\ and the "radial mode contribution" to 
the Schwinger-Dyson equation is not infrared divergent. 
We note, however, that the radial mode contribution is 
subleading as compared to the Goldstone mode contribu- 
tion such that our approximation catches the dominant 
features for the behavior of A^. 

Finally, in the superfluid phase the potential (jElip also 
contains a cubic term oc 0i02 ^ ^fh!^ ■ This term con- 

2 

tributes to the Schwinger-Dyson equation for vrv^ . The 

coefficient of this contribution oc A^ A^po vanishes for 
A^ = 0. Nevertheless, for a momentum dependent A^ 
one has to take it into account, as well as similar cor- 
rections to the Schwinger-Dyson equation for A^ which 
involve quintic couplings. 

APPENDIX F: NUMERICAL PROCEDURES 

In this section we give a short summary which equa- 
tions are actually used for numerical solutions. All quan- 
tities are given in dimensionless renormalized units. We 
first give a complete list relating dimensionful, dimen- 
sionless and dimensionless renormalized parameters, cou- 
plings and fields. 

(i) Relations between dimensionful and dimensionless pa- 
rameters (the Fermi energy is given by = A;|,/2Af) 

q — qjk, T^T/ep, v — v/ep, c = akF, 

a = a/ep, = 2Mh^/k^p'^ , 

A^ = 2MA^, = fn1/eF, A^ = 2MkF\^ 

01 = kp^^'^Tp, = kp^^'^cf), 

p ^ 4>*4> = kp^p, f^hlp^r/e%. (Fl) 

(ii) Relations between dimensionless and dimensionless 
renormalized parameters 

^0 = A^/Z,f,, m^^m^/Z^, v = i)/Z^, 

A0 = \4,lZl, h^ = h^/Z^J^, p^Z^p. (F2) 

All other parameters, couplings and fields are invariant 
under a rescaling with Z^. In particular, note the in- 
variance of the concentration c and the dimensionless 
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superfluid gap parameter f. In our approximation, we 
neglect details of the renormalization coefficients for the 
Matsubara frequencies and use 



Z^i^t/Z^ — 1. 



(F3) 



The input parameters are the detuning !>, the tempera- 
ture T and the Yukawa coupling h^, with c following from 
eq. ([M]) . Alternatively, we could choose c~^, T. We 
have to solve the field equations for the effective chemical 
potential a (density equation) and the field expectation 
value 4>. The latter has the general form 



du 
dp 



(po) • (/>o = 



(F4) 



In the superfluid phase it determines the expectation 
value 00- In the normal or symmetric phase the field 
expectation value vanishes, 0o = 0. The field equation 
for 4> is now replaced by the gap equation for the boson 
mass term = du{0)/dp. The two equations provide a 
closed system, whose solution for a and po or deter- 
mines all observables of the crossover problem. 

In the superfluid or symmetry-broken phase we have 
4>o ¥^ Q such that du{po)/dp = must hold. The density 
equation and the latter condition are then solved for a 
and f — h?^(f)*(f). However, the full four-boson couphng 
enters both the density equation and the gap equation. 
Therefore, an additional gap equation for is needed, 
cf. sect. I Villi and app. [E] Hence the solution is given 
by CT, f, A^. In the "zero momentum approximation" we 
find A^ ~ 0. 

All other quantities of interest can then be recon- 
structed from the solution. For example, we can obtain 
the bare noncondensed molecule fraction and the conden- 
sate fraction by 



M 



(F5) 



Here Q.m depends on the full boson mass term m^, and 

m^^^^,A^^^ (eqs. (|I35[ If 36p ) have to be evaluated at 
f = 0. 

The field equation for a is equivalent to the condition 

i = nF + nM (F8) 

which uses the dressed density fractions. Here 



ftp 



(27r)3 



(e^-' + l)"', 



(2^)3 V 



(ii) Superfluid phase. 

Now f and a are determined by the equations 



(F9) 



(FIO) 



AF)2 



n 



F,0 



C 



I. 



(FII) 
(FI2) 



These equations are still coupled to a gap equation for 
A^ (cf. app. [E] for details). Here f^M is given by eq. 
(fIT3)) f»Aj = Qn'^kp^nM), by eq. ^ and 



i2F,0 — — StT 



(271)3^0 



tanh7^-l). (FI3) 



(iii) Phase transition. 

At the critical line T = Tc we have 

solve 



and we 



where hm is given by eqs. pi3p . The dressed noncon- 
densed molecule fraction and condensate fraction are the 
given by 



m 



(F)2 



1 

3^ 



ftp + n 



M 



1 



(FI4) 
(FI5) 



(F6) 



Let us now give the explicit formulae used in the dif- 
ferent regions of the phase diagram, 
(i) Normal phase. 

The full boson mass is determined by the condition 



(F)2 



1 

3^ 



M- 



(F7) 



at f = for a and Tc- Here is evaluated at — 
and reads 
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